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Preface 


In  the  Air  Force,  there  is  a  strong  need  for  the  lightest 
weight  electrical  generator  with  the  highest  possible  efficiency. 
This  brought  about  a  thesis  topic.  In  the  topic  it  was  proposed 
that  two  students  investigate  the  effects  of  rotor  geometry 
on  improving  the  harmonic  performance  of  a  synchronous  gen¬ 
erator.  If  the  harmonic  performance  could  be  improved  enough, 
it  was  hoped  that  wave  shaping  circuitry  external  to  the 
generator  could  be  eliminated. 

We  chose  this  thesis  topic  because  of  the  depth  of  our 
early  graduate  courses  in  power  systems  and  our  keen  desire 
to  increase  our  knowledge  in  this  area. 

We  would  like  to  thank  Dr.  Frederick  Brockhurst  for  his 
guidance  throughout  all  stages  of  this  thesis  effort.  We 
would  also  like  to  thank  Mr.  James  Kummeth  for  his  assistance 
on  the  Finite  Element  analysis  computer  work. 

Cecil  Stuerke 
David  W.  Nordquist 
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Abstract 

The  use  of  rotor  geometry  as  a  specification  in  the 
control  of  generator  voltage  and  flux  density  harmonics  is 
studied.  The  possibility  of  generating  a  specific  voltage 
waveform  at  the  load  with  only  a  three-phase  rectifier  between 
the  AC  generator  and  the  load  is  investigated. 

Since  the  specified  load  voltage  is  DC,  the  internally 
generated  phase  voltage  waveform  must  be  one  which  keeps  con¬ 
stant  the  greatest  instantaneous  potential  difference  between 
any  two  phases. 

The  armature  reaction  NIMF  is  expressed  as  an  infinite  set 
of  traveling  waves  (with  respect  to  the  stator)  with  an  infinite 
set  of  velocities.  It  is  shown  that  the  rotor  NIMF  cannot  can¬ 
cel  armature  reaction  everywhere  because  the  rotor  NIMF  has 
only  the  velocity  of  the  rotor  (not  an  infinite  set  of  velo¬ 
cities)  . 

The  MMFs  present  at  all  the  stator  slots  of  an  armature 
coil  group  are  reflected  into  a  single  slot  in  an  attempt  to 
control  the  air  gap  flux  density  harmonics  by  selectively 
positioning  increments  of  rotor  reluctance.  This  method 
fails  because  a  general  solution  for  the  reflected  flux 
densities  (appropriately  synchronized  with  each  other)  could 
not  be  obtained.  This  method  can  be  further  investigated  if 
such  a  general  solution  can  be  found. 
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Specifying  the  air  gap  reluctance  incrementally  as  the 
rotor  moves  under  successive  stator  slots  is  also  tried  as  a 
means  of  harmonic  control.  This  method  proves  unworkable 
due  to  negative  NIMF  requirements  and  conflicting  geometry 
specifications . 

Finite  Element  analysis  is  intioduced  as  a  computer  aid 
to  the  design  of  a  magnetic  structure.  A  specific  computer 
program,  MAGNETIC,  is  introduced.  Another  computer  program 
is  developed  which  generates  input  data  for  MAGNETIC.  This 
input  data  is  for  the  pole  shoe  and  air  gap  regions  of  the 
generator.  The  program  allows  for  a  variable  air  gap  shape. 
Recommendations  are  included  which  may  allow  the 
successful  specification  of  rotor  geometry  as  a  means  of 
harmonic  control,  or  which  may  prove  the  concept  impossible. 
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THE  EFFECTS  OF  ROTOR  GEOMETRY 


ON  THE  HARMONIC  PERFORMANCE  OF 
SYNCHRONOUS  GENERATORS 

I.  Introduction 


Background 

For  more  than  a  decade,  synchronous  AC  generators  have 
been  used  to  supply  power  to  loads  that  require  other  than 
a  pure  sinusoidal  AC  voltage.  The  circuit  configuration  of 
a  typical  example  is  shown  in  Figure  1. 

The  AC  generator,  also  called  an  alternator,  allows 
sliprings  to  replace  troublesome  DC  generator  commutators. 

In  the  case  of  brushless  AC  generators,  sliprings  are  even 
eliminated.  In  high  voltage  applications,  AC  generators 
have  the  additional  advantages  of  a  more  easily  insulated 
stationary  armature  (Ref  6:336). 

In  order  to  supply  a  non-sinusoidal  voltage,  a  conven¬ 
tional  AC  generator  requires  external  wave  shaping  circuitry. 

A  DC  load  voltage  requires  that  the  AC  generator  output  be 
rectified  and  also  requires  a  large  filter.  A  DC  generator 
application  may  use  a  smaller  filter,  but  this  does  not 
outweigh  the  commutator  disadvantage. 

The  relative  simplicity,  reliability,  and  insensitivity 
to  wear,  corrosion,  and  breakage  make  AC  generators  beneficial 
even  in  DC  applications  (Ref  2:76)  and  (Ref  7).  An  AC  gener¬ 
ator  with  a  tailored  voltage  output  adds  to  these  advantages. 
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The  voltage  in  each  AC  generator  coil  varies  in  time, 
but  is  not  automatically  sinusoidal.  The  voltage  wave  shape 
depends  on  the  magnetic  flux  density  field  in  the  generator's 
air  gap.  The  field  is  distorted  due  to  the  effects  on  salient 
pole  geometry,  conductor  slots,  and  armature  reaction.  The 
time  variation  of  the  distorted  voltage  wave  shape  may  be 
resolved  into  the  sum  of  fundamental  and  higher  order  Fourier 
components . 

The  harmonics  of  the  voltage  time  variation  must  be 
controlled.  The  first  few  harmonics  contain  significant 
power  that  the  prime  mover  must  supply.  If  the  power  from 
some  of  these  harmonics  must  be  trimmed  from  the  load,  it  is 
generally  dissipated  as  wasted  heat  in  a  conventional  gener¬ 
ator  and  filter  circuit.  Without  proper  control,  the  voltage 
wave  shape  inside  the  generator  may  have  peaks  that  require 
increased  armature  insulation.  An  undesirably  large  and 
heavy  generator  is  required  to  solve  the  preceding  problems. 

Conventional  generator  design  includes  methods  for  har¬ 
monic  control.  Generator  geometry  is  normally  used  to  sup¬ 
press  the  most  significant  undesired  harmonics.  Specification 
of  the  distribution,  pitch,  and  skew  of  the  armature  coils 
as  well  as  the  rotor  pole  arc  provides  a  direct  and  effective 
means  of  harmonic  control.  If  the  wave  shape  can  be  trimmed 
to  its  fundamental  component,  the  analysis  necessary  to 
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shape  the  waveform  to  the  needs  of  the  load  is  much  easier. 

For  any  waveform,  the  shaping  amounts  to  the  reintroduction 
of  properly  sized  higher  harmonics  at  the  expense  of  the 
fundamental. 

When  the  generator  is  dedicated  to  specific  load  which 
requires  other  than  the  fundamental  frequency  voltage  and 
current,  the  conventional  design  procedure  is  a  disservice. 

If  the  voltage  wave  shape  is  obtained  directly  at  the  gener¬ 
ator  terminals,  it  is  possible  to  meet  the  load  power  require¬ 
ments  with  a  smaller  filter  and  at  the  same  time  reduce 
generator  losses.  Rotor  geometry  control  can  help  provide 
the  correct  voltage  at  the  generator  terminals. 

Problem  and  Scope 

The  purpose  of  this  thesis  effort  is  to  investigate  the 
possibility  of  generating  a  specific  DC  voltage  waveform  at 
the  load  with  only  a  rectifier  between  the  AC  generator  and 
the  load.  Rotor  geometry  modifications  are  to  be  investigated 
as  the  main  means  of  control.  This  control  method  is  applied 
to  the  generation  of  specific  phase  voltages  (by  controlling 
generator  flux  density  harmonics)  that  will  yield  the  desired 
load  voltage  using  only  rectification  to  modify  the  voltage 
wave  shape.  This  effort  reviews  the  theory  by  which  rotor 
geometry  and  the  voltage  wave  shape  are  related.  From  this 
basis,  processes  for  determining  a  satisfactory  rotor  geometry 
are  studied. 
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Assumptions 

The  attempts  at  rotor  geometry  determination  through 
required  reluctance  do  not  account  for  the  effects  of  finite 
permeability  or  flux  fringing  in  the  generator  magnetic 
structure.  These  effects  arc  considered  negligible,  but  are 
discussed  as  they  are  observed  in  the  finite  element  analysis 
(Section  X).  The  basic  generator  design  and  the  load  para¬ 
meters  used  in  the  development  of  numerical  examples  arc 
assumptions  taken  from  the  generator  design  covered  in 
Appendix  I. 

Approach  and  Presentation 

Section  I  is  an  introduction.  Section  II  reviews  the 
theory  tlirough  whicli  the  generator  magnetic  flux  density 
is  related  to  rotor  geometry.  This  provides  the  basis  for 
harmonic  control. 

In  Section  111,  the  load  voltage  requirements  are  used 
to  select  the  waveform  of  the  phase  voltage  and  current.  The 
spec  if  icat ion  of  the  phase  voltage  waveform  determines  the 
required  magnetic  flux  density  versus  time  in  the  generator. 
The  specification  of  the  phase  current  allows  calculation  of 
the  armature  reaction.  Armature  reaction  MMT  (magnetomotive 
force) is  calculated  in  Section  IV. 

Section  V  develoiis  an  expression  for  tlic  rotor  MMl . 
Section  VI  combines  tiic  armature  reaction  and  rotor  MMl'  in 
an  expression  for  the  MMl-  at  the  slots  of  a  phase  group  ver¬ 
sus  time.  An  cxaiiqile  is  included. 
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Section  VII  discusses  the  limitations  in  the  use  of 
rotor  MMF  to  cancel  certain  harmonics  of  armature  reaction. 

Section  VIII  describes  an  attempt  to  control  the  flux 
density  harmonics  with  the  introduction  of  selectively  posi¬ 
tioned  increments  of  rotor  reluctance.  Reasons  for  the 
inapplicability  of  this  approach  are  discussed. 

Section  IX  describes  the  incremental  specification  of 
rotor  reluctance  as  it  is  rotated  into  position  under  an 
armature  coil.  The  insufficiency  of  this  approach  is  also 
discussed. 

Section  X  looks  at  the  effects  of  finite  permeability 
and  flux  fringing.  Finite  Element  analysis  is  shown  to  be  an 
excellent  tool  in  fine  tuning  generator  design  to  account 
for  these  effects. 

Section  XI  recommends  other  approaches  to  the  solution 
of  the  harmonic  control  problem  for  further  research. 
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I I .  Theory 

This  section  reviews  the  theory  through  which  the 
generator  magnetic  flux  density  is  related  to  rotor  geometry. 
This  discussion  provides  a  basis  for  the  control  of  the  flux 
density  harmonics. 

In  electromechanical  energy  conversion,  the  magnetic 
field  is  the  intermediate  stage.  A  rotating  generator  shaft 
causes  magnetic  lines  of  flux  generated  by  the  rotor  coils 
to  cut  the  armature  coils.  A  voltage  is  induced  in  each 
coil  according  to  Faraday's  Law 


(/ 


coil 


-d<t>  /dt 


(1) 


where 


V  is  the  instantaneous  induced  voltage 
coil  ® 

<{i  is  the  magnetic  flux  in  webers 

t  is  the  time  in  seconds 


The  coil  voltages  combine  as  coil  group  voltage  which 
combine  as  phase  in  voltages.  The  phase  voltages  are  recti¬ 
fied  to  provide  the  load  voltage. 

"At  a  particular  rotating  speed,  the  instan¬ 
taneous  volts  per  conductor  are  proportional  to 
the  air  gap  flux  density  at  the  conductor.  The 
wave  shape  of  the  conductor  voltage  versus  time 
is  therefore  the  same  as  that  of  the  air  gap  flux 
density  vs  distance  around  the  periphery."  (Ref  3; 

Chap  6,6) 

The  phase  voltage  is  shaped  by  both  the  flux  density  and  the 
layout  connection  of  the  armature  coils. 
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The  shape  o£  the  flux  density  versus  distance  may  be 
described  as  the  sum  of  harmonic  components.  Each  component 
contributes  to  the  phase  voltage.  The  stator  coil's  pitch, 
distribution,  and  skew  determine  multiplier  factors  which 
modify  the  voltage  (Ref  5:44-55). 

When  the  coil  width  is  less  than  the  width  of  a  rotor 
pole,  its  voltage  is  modified  by  a  pitch  factor. 


where 


sin 


V1T  W  . 

T-  7^) 


(2) 


K  is  the  vth  harmonic  pitch  factor 
pv 

V  is  the  space  harmonic  number  being  considered 

W  is  the  winding  pitch  in  radians  or  slots 

T  is  the  first  harmonic  rotor  pole  pitch  in 
P  radians  or  slots 

When  the  coils  are  not  concentrically  placed,  the 
voltage  is  modified  by  a  distribution  factor. 


where 


sin(vti/m) 
q  sin(vTr/qm) 


(3) 


Kj  is  the  vth  harmonic  distribution  factor 
dv 

q  is  the  width  of  a  phase  belt  in  slots 
m  is  the  number  of  internally  generated  phases 
When  the  coils  are  skewed  with  respect  to  the  rotor 
axis,  the  voltage  is  modified  by  a  skew  factor. 


8 


both  i^oiirccs  of  MMl'  rotatinj;  with  respect  to  tlic  stator, 
rlie  field  produced  i^y  tiicsc  MMf-  sources  may  be  controlled 
tlirough  variation  of  tlie  magnetic  circuit  reluctance. 


Figure  2  One  Pole  o£  the  Rotor  and  the  Stator 
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({I  =  MMF/R 


(5) 


where 

NIMF  is  the  magnetomotive  force  in  ampere  turns 
R  is  the  reluctance  in  henries  ^ 

Reluctance  of  the  magnetic  path  is; 

R  =  (6) 

where 

p  is  the  permeability  of  the  material  in  henries/ 
meter 

is  the  path  length  in  meters 

2 

a  is  the  area  in  meters 

Because  the  permeability  for  air  is  much  lower  than  the 
permeability  of  the  generator  steel,  the  air  gap  usually 
dominates  the  reluctance  and  the  NIMF  drop  around  the  magnet¬ 
ic  circuit.  The  air  gap  reluctance  (above  the  rotor)  travels 
with  respect  to  the  stator  as  does  the  MMF  which  forces  the 
flux  through  it. 

If  both  sides  of  Equation  (5)  are  divided  by  the  area  of 
the  flux  tube  (Ref  1:245),  in  which  the  flux  is  found,  it  is 
evident  that  the  flux  density  is  related  to  the  MMF  through 
the  length  of  the  flux  tube. 

B  =  u  MMF/£  (7) 

m 

where 
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2 

B  is  the  magnetic  flux  density  in  webers/meter 
When  the  reluctance  is  due  to  the  air  gap  distance,  this 
becomes 


B  =  a  MMF/6(x) 


(8) 


where 


6(x)  is  the  air  gap  distance  in  meters 

X  is  the  distance  along  the  circumference  of  the 
inside  of  the  stator  in  meters  (stator 
coordinates) 

The  chosen  method  for  tailoring  the  magnetic  flux  density 
is  first  to  calculate  the  MMF  variation  at  each  slot,  and  then 
to  adjust  the  reluctance  through  variation  of  the  length  of 
the  air  gap  at  the  slots.  The  magnetic  flux  density  in  turn 
produces  a  phase  voltage  that  becomes  the  desired  voltage  wave 
shape  when  rectified. 
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III.  Phase  Voltage  and  Current  Wave  Shapes 


This  section  presents  the  theory  and  the  selection  of  an 
acceptable  phase  voltage  wave  shape.  The  current  wave  shape 
which  follows  from  the  voltage  is  also  presented.  The  load 
voltage  requirement  is  the  basic  criteria  on  wliich  the  selec¬ 
tions  are  based. 

A  number  of  different  phase  voltage  wave  shapes  can  pro¬ 
duce  the  desired  load  voltage  wave  shape  after  being  rectified 
Thus,  the  harmonic  content  of  the  generator  magnetic  field 
that  produces  the  voltage  and  current  wave  shapes  is  a  design 
opt  ion . 


Voltage 

The  generator  must  supply  the  required  load  voltage  in¬ 
spite  of  the  generator  transmission  line  impedances.  If,  as 
is  shown  in  Figure  3,  the  load  voltage  is  specified  and  the 
impedances  are  known,  the  required  line  to  line  voltage  wave 
shape  may  be  determined  by  circuit  analysis. 


where 


1/ 

g 


(Z, 


Z^) 


(9) 


1/ 

g 


is  the  internally  generated  line  to  line  voltage 
is  the  required  load  voltage 
is  the  load  impedance  in  ohms 
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I'iyurc  5  PIktsc  (Circuit  Impedance 


is  tiie  internal  generator  impedance  in  ohms 


'T 


is  the  transmission  line  impedance  in  olims 


This  analysis  is  to  be  accomplished  harmonic  by  harmonic  from 
the  components  oT  the  load  voltaj^e.  Remember  tliat  tlie  im|ie- 
dances  are  functions  of  frequency.  It  is  the  phase  voltage 
leading  to  the  composite  that  must  be  induced  in  the 
armature  coils  by  tlie  aggregate  magnetic  field. 

Tlie  shape  of  the  desired  load  voltage  wave  sliape  speci¬ 
fies  the  phase  voltage  difference  (between  any  two  e.xternal 
generator  phases)  during  the  period  in  which  it  is  applied  to 
the  load.  This  is  because  the  three-phase  rectifier  ajiplies 
the  greatest  instantaneous  potential  difference  between  any 
two  piiases  to  tiio  load.  If  the  load  voltage  is  Dt!,  tlie  ma-ximum 
phase  difference  is  required  to  be  constant.  As  tlie  voltage 
in  one  phase  falls  below  the  necessary  value,  another  phase 
must  rise  to  replace  it. 

hqually  loaded  symmetrical  phases  will  have  identical 
voltage  wave  shapes.  Successive  phases  will  be  evenly  dis¬ 
placed  through  2tt  electrical  radians  of  rotor  rotation. 

T.xample  phase  voltage  wave  shapes  which  meet  the  I'C  load 
reepj i rement s  are  shown  in  l  igure  4.  Notice  that  in  figure  4c, 
overlapping  voltages  are  possible.  In  ,  e  cases,  current  is 
supplied  from  more  than  two  phases  at  a  time. 

The  shape  of  the  phase  vo 1 1 ages  as  they  transition  is 
not  specified  by  the  load  voltage  requirements.  It  is 
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unnecessary  for  the  phase  voltage  transition  to  encourage 
proper  current  commutation  because  the  current  in  the  load 
is  controlled  by  the  applied  voltage.  The  rate  of  change 
of  the  phase  voltage  may  be  limited  in  order  to  protect  the 
silicon  controlled  rectifiers  and/or  diodes  used  in  the 
three-phase  rectifier. 

Flexibility  in  the  choice  of  the  phase  voltage  wave 
shape  is  useful  because  this  wave  shape  must  be  controlled 
by  a  properly  engineered  magnetic  flux  density  in  the  air 
gap. 

This  flexibility  may  be  exploited  under  rotor  geometry 
control  in  order  to  accomplish  flux  density  control  with  a 
minimum  mean  air  gap.  This  will  decrease  leakage  and  improve 
generator  performance.  Because  the  phase  voltage  determines 
phase  current,  its  shape  also  affects  the  magnetic  field 
shape  resulting  from  the  armature  reaction. 

The  shape  of  the  load  voltage  during  its  application  to 
the  load  and  its  transition  affects  its  harmonic  content. 

The  proper  choice  of  the  phase  voltage  wave  shape  as  it  is 
applied  to  the  three-phase  rectifier  may  improve  generator 
performance  and/or  simplify  calculations.  Once  the  effective 
rotor  length,  the  rotor  tip  velocity,  and  the  phase  voltage 
are  known,  the  required  magnetic  flux  density  is  also  known. 


coil 


Bvf 


(10) 
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where 


l^coii  is  the  required  coil  group  voltage 

V  is  the  rotor  tip  velocity  in  meters/second 

I  is  the  effective  rotor  length  in  meters 


Voltage  Example 

See  Appendix  I  (Generator  Design)  for  a  reference  on  the 
generator  specifications  used  in  this  example. 

For  the  example  carried  through  this  thesis  effort,  the 
desired  load  voltage  is  constant  (DC) .  Although  the  internally 
generated  phase  voltage  is  properly  derived  from  the  desired 
load  voltage  according  to  Equation  (9) ,  this  example  begins 
with  the  specification  of  the  line  to  line  voltage.  The 
value  of  the  line  to  line  voltage  is  chosen  as  796.74  volts. 
This  is  /T  times  the  RMS  phase  voltage  of  the  sample  sinusoidal 
AC  generator.  It  is  equal  to  the  RMS  voltage  the  generator 
could  supply  to  a  load  if  it  were  wye  connected  and  if  2^  and 

Z  were  negligible. 

R 

The  generator  phase  voltage  output  comes  from  the  series 
connection  of  four  coil  groups  (two  coil  groups  per  pole  pair 
from  the  generator's  two  pole  pairs).  The  selected  internal 
phase  voltage  wave  shape  is  shown  in  Figure  5.  The  rectifi¬ 
cation  applies  the  maximum  potential  difference  between  any 
two  external  phases  to  the  load.  Considering  that  this 
difference  is  across  two  external  phases  and  considering  the 
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four  coil  groups  in  series  per  external  phase,  the  required 
coil  group  voltage  is: 

l'*coii  =  796.74/8  =  99.59  volts  (11) 

where 

A 

1/  -1  is  the  peak  value  of  1/ 

coil  ^  coil 

The  rotor  tip  speed  (v)  is  119.7  meters/second  and  the 
effective  rotor  length  is  .6438  meters.  The  total  flux  den¬ 
sity  which  the  five  conductors  in  each  coil  group  must  cut 
has  the  same  shape  versus  time  as  From  Equation  (10), 

the  maximum  flux  density  is: 

* 

B  =  — =  1.292  webers/meter^  (12) 

max  VC  ^ 

where 

B  is  the  maximum  flux  density  in  webers /meter “ 
max  ^ 

Figure  5  gives  the  voltage  versus  time  (external  to  the 
generator)  for  the  three  phases.  Phase  B  lags  phase  A  inter¬ 
nally  in  the  generator  by  Tr/3  electrical  radians  (60  electri¬ 
cal  degrees).  Because  phase  B  is  connected  in  reverse,  it 
lags  phase  A  by  an  additional  ti  electrical  radians  (180 
electrical  degrees)  for  a  total  lag  of  4tt/3  electrical 
radians  (240  electrical  degrees).  The  flux  density  for  each 
phase  has  the  same  shape  as  the  phase  voltage  variation. 

By  Fourier  analysis,  (Ref  14:Chap  2)  the  harmonic  com¬ 
ponents  of  the  desired  total  flux  density  per  phase  are  given  by 
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Figure  5  Time  Variation  of  the  Chosen  Internal  Phase 
Voltage  and  Flux  Density 


B  =  -  I 
^  g=l 


®max  g  sin  (gut) 


(13) 


where 

u  is  the  frequency  in  radians/second 
g  is  the  time  harmonic  number  being  considered 


Current 

The  phase  current  versus  time  is  determined  by  the  shape 
of  the  phase  voltage  waveform  as  well  as  the  circuit  impe¬ 
dances.  Between  current  commutation  periods,  the  phase 
current  is  the  specified  load  voltage  divided  by  the  speci¬ 
fied  load  impedance.  The  value  of  the  phase  current  during 
commutation  may  be  determined  by  the  use  of  harmonic  compon¬ 
ent  analysis  (Ref  12:400).  Because  the  circuit  may  have 
frequency  dependant  impedances,  each  harmonic  component  of 
the  current  may  have  a  different  angle  by  which  it  lags  the 
voltage  of  the  same  harmonic.  The  composite  waveform  accounts 
for  the  current  lag.  Another  approach  is  to  analyze  the  com¬ 
posite  voltage  waveform  directly  (Ref  7). 

Current  fxample 

See  Appendix  I  (Generator  Design)  for  a  reference  on  the 
generator  specifications  used  in  this  example. 

Because  it  is  necessary  to  maintain  the  load  power 
at  100  kilowatts,  the  DC  load  current  is  125.51  amperes  (100 
kilowatts  divided  by  460*/T  volts).  There  are  no  parallel 
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current  paths  inside  the  generator.  Therefore,  between 
commutation  periods,  the  phase  current  is  also  125.51  amperes. 

Figure  6  shows  the  circuit  of  Figure  1  during  a  phase  B 
to  a  phase  A  commutation.  Phase  C  is  conducting  throughout 
this  time,  so  no  diodes  are  shown.  It  is  assumed  that  the 
transmission  line  impedance  and  the  generator  resistance  are 
negligible. 

The  phase  voltage  wave  shape  was  given  in  Figure  5.  A 
constant  voltage  (4  •  applied  across  the  inductive 

coils  of  phases  A  and  B.  The  current  in  phase  B  linearly 
increases  from  zero  and  the  current  in  A  linearly  decreases 
toward  zero.  The  generator  phase  inductance  for  this  example 
is  chosen  such  that  the  current  commutation  angle  is  tt/9 
electrical  radians  or  equivalently  20  electrical  degrees. 

The  current  wave  shape  of  the  three  phases  is  given  in 
Figure  7.  Notice  how  they  follow  the  voltages  of  Figure  5. 

The  time  variation  of  the  current  for  the  three  phases  at 
selected  times  is  given  in  Table  1. 

The  harmonic  components  of  the  phase  current  are: 

I(t)  =  I  (36  il  [cos(gn/9)  -  cos(g27T/9)] 

g=l 

+  i  [sin(g2iT/9)  -  sin(gTi/9)]}  sin  g(wt-iT/18)  (14) 

where 

I(t)  is  the  phase  current  in  amperes 

I  is  the  maximum  value  of  the  phase  current 

max 
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Table  1  Normalized  Current  in  Six  Internal 
Phases  versus  Time 


TIME 

PHASE 

(wt) 

A 

B 

C 

A’ 

B' 

C’ 

-  Stt  / 1 5 

-1 

1 

0 

1 

-1 

0 

-47r/15 

-1 

1 

0 

i. 

- 1 

0 

-3tt/15 

-1 

1 

0 

1 

-1 

0 

-Ztt/IS 

-7 

•  / 

1 

-  .3 

.  7 

-1 

.  3 

-tt/15 

- .  1 

1 

-.9 

,  1 

- 1 

.  9 

0 

0 

1 

-1 

0 

-1 

1 

r/15 

0 

1 

-1 

0 

-1 

1 

2tt/15 

0 

1 

-1 

0 

-1 

1 

3Tr/15 

.  3 

.  7 

-1 

-  .3 

-  .  7 

1 

4ir/15 

.9 

.  1 

-1 

-  .9 

-  .  1 

1 

5tt/15 

1 

0 

-1 

-1 

0 

1 

6tt  / 1 5 

1 

0 

-1 

- 1 

0 

1 

777/15 

1 

0 

-1 

-1 

0 

1 

Stt/IS 

1 

-.3 

-  .  7 

-1 

.  3 

.  7 

97i/15 

1 

-.9 

-  .  1 

-1 

.  9 

.  1 

IOtt/IS 

1 

-1 

0 

-1 

1 

0 

IV.  Armature  Reaction 


Tfiis  section  develops  an  expression  for  the  time  varia¬ 
tion  of  the  MMF  at  each  of  the  stator  slots  due  to  armature 
react  ion. 

As  magnetic  lines  of  flux  inside  a  generator  cut  armature 
coils,  a  voltage  is  induced  at  the  armature  terminals.  When 
a  load  is  connected  to  the  generator,  current  will  flow 
through  the  armature  coils  generating  an  MMF  known  as  armature 
reaction.  Armature  reaction  MMF  opposes  the  MMF  due  to  the 
rotor  field. 

Iv'hen  the  MMF  from  the  rotor  coils  is  added  to  the  armature 
reaction  MMF,  the  total  MMF  is  found.  The  desired  magnetic 
flux  density  can  be  produced  provided  that  the  correct  air 
gap  reluctance  is  present. 

Single- Pliase  Winding 

A  single  coil  from  each  pole  of  a  single-phase  winding 
produces  a  square  MMF  wave  around  the  air  gap  circumference. 

If  the  coils  span  the  full  width  of  a  pole  (Tp)  ,  each  coil 
has  a  width  of  rr  electrical  radians.  The  magnitude  of  the 
MMF  wave  will  oscillate  with  the  frequency  of  the  current  in 
the  coils. 

The  MMF  under  the  armature  coils  may  be  expressed  as 
the  product  of  the  time  variation  of  the  coil  current  and 
the  Fourier  components  of  the  space  variation  of  the  square 
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MMF  envelope.  If  the  x  coordinate  (the  distance  along  the 
circumference  of  the  inside  of  the  stator)  is  zero  at  a  posi¬ 
tion  which  corresponds  to  the  center  of  one  of  the  full 
pitch  coils,  the  space  variation  of  the  MMF  per  conductor  is 
given  by: 


MMF (for  one  conductor) =  I 


max  IT  ^  1  V 
V  =  1 


1 

^  -sin(vTT/2)  cos(vtTx/T  )  (15) 

_  1  V  p 


Because  of  the  term  sin(vTT/2),  the  spacial  variation  has 
only  odd  harmonics.  Because  of  the  symmetry,  MMF(-x)  =  MMF(x), 
the  coefficient  of  the  sin(viTx/Tp)  terms  in  Equation  (15)  are 
all  zero  (and  not  shown). 

The  current  in  the  coils  is  the  phase  current  as 
described  in  Section  III. 


I(t)  =  I 


I  (a  sin  gut  +  b  cos  gut) 


where 

a  are  the  coefficients  of  the  sine  current  harmonics 
g 

b  are  the  coefficients  of  the  cosine  current  harmonics 
g 

Because:  a  sin(ut  +  <t))  +  b  cos(ut  +  <j))  =  a'  sin  ut  +  b'  cos  ut 

a  proper  set  of  coefficients  a  and  b  may  be  established 

g  g 

that  account  for  the  current  lag  angle  of  each  harmonic. 

If,  as  is  usually  the  case,  the  negative  portion  of  the  period¬ 
ic  current  wave  is  a  mirror  image  of  the  positive  position, 
only  odd  harmonics  exist  in  the  time  variation  of  the  current. 


A  generator  phase  may  contain  several  conductors  dis¬ 
tributed  over  one  pole  pitch  in  a  single  coil  group.  Since 
tlie  current  in  each  conductor  is  identical,  the  resultant 

MMl-  mav  be  summed  using  t  he  d  i  s  t  r  i  but  i  on  f  ac  t  or  K,  .  The 

dv 

space  variation  then  becomes: 


MMF(x) 


MMF 

max 


2 

TT 


00 


I 

v  =  l 


-  sin(vTT/2) 


K,  COSfvTTX/?  ) 
dv  ^  '  p-' 


(17) 


where 


NWF  is  the  maximum  value  of  the  NIMF 

max 


MMF 


max 


qs  I 


max 


(18) 


where 

s  is  the  number  of  turns  of  a  coil  group  in  a 
single  stator  slot 

Notice  that  the  coordinate  x  is  now  taken  to  be  zero  in  the 
center  of  a  coil  group. 

The  space  and  time  variation  of  a  single  phase  is  the 
product  of  the  individual  variations. 


MMF(x,t) 


MMF 


max 


i  I  [a  sin(gu)t)  +  b  cos(ga)t)] 


‘  ^  I  ^  sin(vn/2)  COS(vTrx/T  ) 

"  v  =  l  ^ 


(19) 


where 

•  indicates  a  product 

g  is  the  time  harmonic  number  being  considered 
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Distribution  lixample 

See  Appendix  I  (Generator  Design)  for  a  reference  on 
the  generator  specifications  used  in  this  example. 

The  distribution  factor  for  this  example  is: 


s in  (uTi  /m)  _  sin(v^/b) 

q  sin(yTr/mq)  5  sin(vTT/30)  ^  ' 

Polyphase  Windings 

If  there  is  more  than  one  coil  group  above  each  pole, 
the  phase  of  each  coil  group  must  be  considered  in  calculating 
the  resultant  MMP.  The  current  in  each  coil  group  is  depend¬ 
ent  upon  the  connection  of  the  internal  machine  phases  to  the 
load.  If  there  arc  m  internal  phases,  there  will  be  m  possible 
different  time  variations  of  current  in  the  stator.  Hach 
current  wave  shape  will  be  equally  distributed  in  time  with 
a  separation  angle  of  27t/m  electrical  radians. 

As  the  current  reaches  maximum  in  successive  phases, 
the  maximum  fundamental  MMP  occurs  at  the  center  of  the  suc¬ 
cessive  coil  groups.  Successive  coil  groups  correspond  to 
successive  internal  machine  phases.  The  M.\I1  s  are  therefore 
distributed  Zn/m  electrical  radians  apart. 

The  MMF  due  to  all  of  the  phases  combined  is  the  maximum 
NIMF  times  the  product  of  the  time  K(t,r)  and  space  K(x,r) 
variations  for  all  of  the  m  internal  phases. 


K(t,r)  =  'l  [a  sin  g(ut  +  rZir/m)  +  b^  cos  g(a)t  +  rlTr/m)]  (211 
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where 


K(t,r)  is  the  time  MMF  variation 
r  is  the  individual  internal  machine  phase 

oo 

"  7  I  t  ^  sin(vTr/2)  K,^  cos  v  +  r2n/m)]  (22) 

\;  =  1  p 


where 

K(x,r)  is  the  space  MMF  variation 

m- 1 

MMF(x,t)  =  MMF^^^  I  [K(t,r)  •  K(x,r)]  (23) 

r=0 


Example 

See  Appendix  I  (Generator  Design)  for  a  reference  on  the 
generator  specifications  used  in  this  example. 

The  space  variation  term  in  Equation  (22)  converges  as  it 
is  summed  over  the  integer  v  from  one  to  infinity.  The 
resulting  value  of  the  expression  for  the  six  internal  phases 
(r  =  0  to  r  =  5)  and  for  the  position  of  the  five  slots  in 
phase  C  (Figure  2)  is  given  in  Table  2. 

Table  2  Normalized  Spacial  Component  of 

Armature  Reaction,  K(x,r),  at  Slots 
of  Phase  C  (Figure  2) 


Phase 

B 

C 

B' 

B 

BOB 

0 

1 

2 

3 

4 

5 

Tr/30 

3ti/30 

5tt/30 

7n/30 

97t/30 

.5  .5  -.4  -.5  -.5  .4 

.5  .5  -.2  -.5  -.5  .2 

.5  .5  0  -.5  -.5  0 

.5  .5  .2  -.5  -.5  -.2 

.5  .5  .4  -.5  -.5  -.4 
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IVith  the  current  wave  shape  that  was  established  in 
Section  III,  it  is  possible  to  evaluate  the  armature 
reaction  MMF  at  each  of  the  slots  versus  time.  According  to 
liquation  (23)  and  using  data  from  Table  1  and  Table  2,  the 
products  of  the  time  variation  and  the  space  variation  may 
be  summed  over  the  applicable  phases.  The  number  of  internal 
phases,  m,  is  equal  to  six  in  this  example. 

■Notice  that  both  the  time  and  the  space  variation  terms 
for  phase  A'  are  opposite  in  sign  from  the  entries  for  phase 
A.  The  product  of  the  space  and  time  variation  terms  leads 
to  a  magnetic  field  from  A'  which  is  equal  to  the  magnetic 
field  from  A  both  in  magnitude  and  sign.  Likewise  the  phases 
B  and  B'  and  the  phases  C  and  C  produce  similiar  results. 

From  Figure  2,  a  positive  current  in  phase  B  produces  an  MMF 
under  the  slots  of  phase  C  which  is  opposite  in  direction  to 
the  MMF  produced  by  a  positive  current  in  phase  A. 

The  entries  in  Table  3  (Total  Normalized  .Armature  Reaction 
MMF)  were  obtained  by  changing  the  sign  of  the  phase  R  and  B' 
currents  in  Table  1  and  by  multiplying  the  resultant  Tabic 
1  array  by  the  transpose  of  the  array  in  Table  2.  This 
multiplication  is  the  direct  normalized  interpretation  of 
Equation  (23).  In  order  to  express  the  actual  .MMF,  the 
entries  must  be  multiplied  by  both  the  total  number  of  con¬ 
ductors  in  series  per  internal  machine  phase  and  by  the 
current  which  each  conductor  carries.  Put  another  way,  the 
entries  in  Table  3  multiplied  bv  MMI  (Equation  18)  will 
yield  the  actual  MMF. 
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V.  Rotor  Mf.lF 


This  section  develops  an  expression  for  the  MMF  due  to 
the  rotor  coils. 

The  rotor  coils  are  a  significant  source  of  MMF  in  the 
generator.  Therefore,  the  affect  on  the  controlled  magnetic 
circuit  must  be  accounted  for. 

Two  different  effects  from  the  rotor  are  to  be  investi¬ 
gated  in  the  search  for  a  solution  to  the  problem  of  harmonic 
control.  In  Section  VII,  the  approach  is  to  use  the  rotating 
rotor  field  MMF  (rotating  with  respect  to  the  stator)  to 
cancel  certain  rotating  MMF  (again  rotating  with  respect  to 
the  stator)  components  created  by  the  armature  reaction. 

This  procedure  requires  a  properly  designed  MMF  over  the  rotor 
surface . 

In  contrast.  Section  VIII  and  Section  IX  capitalize  on 
the  reluctance  of  the  rotor  for  control.  For  this  case, 
the  assumption  is  made  that  most  of  the  MMF  drop  is  across 
the  air  gap  reluctance  and  that  the  rotor  MMF  may  be  repre¬ 
sented  (first  order)  as  a  square  wave  over  the  rotor  surface. 
With  these  assumptions,  the  development  of  the  rotor  MMF 
expression  and  its  translation  into  stator  coordinates  is 
cont inued . 

Without  fringing,  the  width  of  each  square  MMF  pulse  is 

equal  to  the  rotor  pole  arc  bp.  Therefore,  the  pitch  factor 

K  ,  may  be  used  to  form  the  MMF  expression, 
pv 
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MMF 

rotor 


— s  It:  I  —  sin  (virb  /2t  )  cos  (vttx'/T  )  (24) 

•FT  r  F  V  P  P  ^ 


where 


MMF  ^  is  the  rotor  MMF  in  ampere  turns 
rotor  ^ 


s  is  the  number  of  turns  in  a  rotor  coil 

r 

Ip  is  the  rotor  current  in  amperes 


bp  is  the  rotor  pole  arc  in  radians  or  slots 

x'  is  the  distance  along  the  outside  circum¬ 

ference  of  the  rotor  in  meters  (rotor 
coordinates) 


The  zero  (origin)  of  the  rotor  coordinate  system  is  taken  at 
a  pole  center. 

The  MMF  expression  for  non-constant  MMF  across  the  pole 
face  may  also  be  achieved  by  Fourier  analysis.  In  this  case, 
saturation  and  fringing  are  the  means  of  achieving  the  desired 
effect.  Finite  Element  analysis  (Section  X)  is  an  important 
tool  to  use  when  considering  saturation  and  fringing. 

Each  harmonic  component  of  the  rotor  MMF  travels  with 
respect  to  the  stator  and  has  the  same  velocity  as  the  rotor. 
The  MMF  due  to  the  rotor  can  be  transformed  into  the  stator's 
set  of  polar  coordinates  by  transforming  each  harmonic  compon¬ 
ent  into  a  forward  traveling  wave. 

00 

MMF  =  -s  It'  y  —  sin  (vitb  /2t  )  cos  [  (vttx/t  j-vajt] 

rotor  7T  r  F  v  ^  p  p"^  ‘  ^  p^  ' 

(25) 
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The  MMF  which  drives  a  flux  tube  may  be  moved  around 
that  tube  without  changing  the  flux  density.  Therefore  with 
radial  flux  lines,  the  traveling  square  wave  can  be  projected 
across  the  air  gap  with  the  same  angular  coordinates. 

As  was  true  of  the  armature  M'lF  expression,  the  rotor 
.^IMF  versus  time  can  be  evaluated  at  any  x  by  substitution  of 
that  value.  Rotor  MMF  is  opposed  by  the  armature  reaction 
MMF.  The  'IMFs  may  be  algebraically  summed  to  obtain  an 
expression  for  the  total  MMF  along  the  stator  (Section  \']  )  . 

F.xample 

See  Appendix  I  (Generator  Design)  for  a  reference  on  the 
generator  specifications  used  in  this  example. 

The  chosen  pole  arc  is  2tt/3  electrical  radians.  There 
are  533  ampere  turns  in  each  rotor  winding.  The  rotor  MMF  is 
a  square  wave  which  travels  with  the  rotor  under  the  armature 
coils.  The  time  variation  of  the  rotor  MMF  under  the  slots 
of  phase  C  (Figure  2)  is  given  in  Table  4. 


Table  4  Rotor  MMF 


T I  Mi;  \  X 
(ait)  \ 

11/30 

3ii/30 

5i:  /30 

771/30 

9it/30 

-57i/15 

0 

0 

0 

0 

0 

-47i/15 

533 

0 

0 

0 

0 

-3n/15 

5  3  3 

533 

0 

0 

0 

-27r/15 

533 

533 

533 

0 

0 

-ti/15 

533 

5  35 

533 

533 

0 

0 

533 

535 

533 

533 

553 

tt/15 

533 

553 

535 

533 

533 

Ztt/IS 

533 

533 

53  5 

535 

55  3 

3V15 

5  35 

5  3  5 

533 

553 

533 

4h/15 

533 

535 

535 

533 

533 

5ti/15 

5  3  3 

533 

5  33 

533 

555 

6tt  / 1  5 

0 

555 

555 

5  33 

53  3 

7ti/15 

0 

0 

555 

5  33 

5  3  3 

8tt/1  5 

0 

0 

0 

533 

5  3  5 

9tt/15 

0 

0 

0 

0 

533 

1  0  n  /  1  5 

0 

0 

0 

0 

0 

VI.  Total  MMH 


Section  IV  developed  an  exjiression  for  the  armature 
reaction  MMl-  and  Section  \’  developed  an  expression  for  the 
rotor  MMF.  These  MMFs  arc  superimposed  in  the  generator's 
air  gap.  F.ach  MMF  produces  a  magnetic  flux  density  with 
space  and  time  variations  that  are  superimposed.  This  section 
develops  an  expression  for  the  total  MMF. 

At  first  glance  it  seems  possible  to  calculate  the  flux 
density  due  to  armature  reaction  and  then  to  adjust  tlic  rotor 
reluctance  so  the  rotor  MMF  produces  a  cancelling  flux  den¬ 
sity.  This  idea  is  flawed  because  a  change  in  the  magnetic 
circuit  also  changes  tlie  magnetic  flux  density  driven  by  the 
armature  reaction  M'lF.  Therefore,  instead  of  the  superposi¬ 
tion  of  flux  density,  it  is  the  superposition  of  the  MMF 
which  must  be  accomplished.  The  desired  flux  density  may 
then  be  obtained  from  the  total  MMF  by  adjusting  the  air  gap 
reluctance . 

The  total  MMF  in  the  generator  is  the  sum  of  the 
expressions  in  liquation  (25)  and  liquation  (25).  Because  of 
the  direction  defined  for  positive  current  (Figure  2),  the 
armature  reaction  MMF  and  the  rotor  MMl-  add. 

As  is  shown  in  .Section  VII,  specification  of  rotor  MMF 
is  not  sufficient  to  control  the  generator  MMl-  for  all  of  the 
air  gap  for  all  times.  The  armature  reaction  MMF  produces 
both  space  and  time  variations  that  are  not  controlled  by  tlic 


rotor  NLMF.  Substitution  of  the  space  or  the  time  into  the 
expressions  however,  brings  the  terms  of  both  the  armature 
reaction  MMF  and  the  rotor  MMF  to  a  countably  infinite  series. 
The  coordinates  of  the  stator  conductors  are  important  since 
it  is  at  these  points  that  the  time  variation  of  the  flux 
density  creates  the  desired  voltage  wave  shape.  hith  this 
substitution,  control  of  the  harmonics  is  approached  in 
Section  VIII  and  Section  IX. 

Example 

See  Appendix  I  (Generator  Design)  for  a  reference  on 
the  generator  specifications  used  in  this  example. 

The  actual  armature  reaction  MMF  is  obtained  by  multi¬ 
plying  the  normalized  armature  reaction  MMF  in  Section  IV 
by  MMF  (Equation  18) .  The  total  MMF  can  be  obtained  by 
summing  the  armature  reaction  MMF  and  the  rotor  MMF.  This  is 
equivalent  to  multiplying  Table  3  by  ^nd  adding  the 

result  to  Table  4.  The  total  MMF  is  shown  in  Table  5. 
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Table  5  Total  MMF  at  Phase  C  Slots  versus  Time 


TImK  X 
(wt)  ^ 

77/30 

377/  30 

577/30 

777  /  3  0 

977/  30 

-Sr/lS 

-  1  2  55.1 

-1255.1 

-1255.1 

-  1255  .  1 

-1255.1 

-  477  /1  5 

-  722.1 

-1255.1 

-  1255  .1 

-  1  2  55.  1 

-1255.1 

-  377  /  1  5 

-  722.1 

-722.1 

-1255 . 1 

-  1255  .  1 

-1255.1 

-  277  /  1  5 

-  385.2 

-458.5 

-533.8 

-1142.1 

-1217.4 

-tt/I  5 

294.5 

68.6 

-157.3 

-333.2 

-1142  .  1 

0 

407  .  5 

156.5 

-94.6 

-  345.6 

-  596.6 

77/15 

407.5 

156.5 

-94.6 

-345.6 

-596.6 

277  /  1  5 

407 , 5 

156.5 

-94.6 

-345.6 

-596.6 

377  /  1  5 

784.0 

533.0 

282.0 

31 . 0 

-220.0 

4it/15 

1537.1 

1286.1 

1035.0 

784 . 0 

5  3  3.0 

577/1  5 

1662.6 

1411.6 

1160.6 

909.5 

658.5 

677/1  5  ' 

1129.6 

1411.6 

1160.6 

909.5 

658 . 5 

777  /  1  5 

1129.6 

878.6 

1160.6 

909.5 

658. 5 

877/  1  5 

1167 . 2 

991.5 

815.8 

1173.1 

997.4 

977/  1  5 

1242.5 

1217.4 

1192.3 

1167.2 

1675.1 

1  077/1  5 

12  55.1 

1255.1 

1255.1 

1255.1 

1  2  55.1 

39 


VII.  MMF  Cancellation 


This  section  considers  rotor  MMF  for  the  cancellation  of 
selected  armature  reaction  MMF  harmonic  components.  For  this 
analysis,  it  is  useful  to  express  both  the  armature  reaction 
MMF  and  the  rotor  MMF  as  traveling  waves. 

Armature  Reaction  MMF 

Equation  (23)  is  tlie  expression  for  the  armature  reaction 
M^iF.  Symmetry  about  the  time  origin  causes  the  cosine  terms 
in  K(t,r)  (Equation  21)  to  be  zero.  The  cosine  term  due  to 
asymmetry  is  subsequently  considered.  For  a  given  g  and  v, 
several  constants  from  Equation  (23)  may  be  combined  by  letting: 

A  =  —  qs  I  a  sin  (v7t/2)  K,  /v  (26) 

gv  7T  ‘  max  g  ^  ^  dv 

wliere 

A  is  a  constant  whose  value  depends  on  g  and  v 

Cl 

Remembering  that  the  cosine  coefficient  in  Equation  (21)  is 
zero,  the  armature  reaction  MMF  expression  from  Equation  (23) 
can  be  rewritten  as: 

CO  00  fil-  1 

MMF(x,t)  =  I  I  I  A  sin  g  (wt +  r2Tr/m) 
g=l  v=l  r=0  * 

•  cos  v[(Trx/Tp)  +  (r2Tr/m)]  (27) 

Using  the  trigonometric  identity 

sin(a)  cos(f)  =  \  sin(a  -  B)  +  ^  sin(nt  +  P)  (28) 


4  0 


Equation  (27)  may  be  rewritten  as: 

OO  QO 

MMr(.,t)  =  V  ;  iA  . 

0=1  ^=1  ® 

|sin  [(o^^t)  +v:Tx/Tp] 

+  sin  [(g'^t)  +  (VTTX/Tp)  -  (g  +  v)^] 

+  sin  [(gujt)  +  (yirx/Tp)  -  (g  +  v)^] 


+  sin 

[  (g'^t) 

-  (VTTX/Tp)] 

2-' 

+  sin 

[  (g(^t) 

-  (v^x/Tp)  -  (g-v)  — 

+  sin 

[  (g<^t) 

-  (VTTX/Tp)  -  (g-v)^ 

+  sin[  (gwt)  +  (vTTx/Tp)  -  (g  +  v)i2^^^i^]+sin[  (gcot)  -  (vrx/Tp)  -  (g-v)  "  "] 

(29) 

Equation  (29)  contains  the  positive  direction  (terms  in  the 
right  column  of  the  equation)  and  negative  direction  (terms 
in  the  left  column  of  the  equation)  traveling  waves. 

When  is  a  multiple  of  m. 


2  TT 

sin  0  +  sin  (0  +  — )  +  sin  (0 

m 


,  4it.. 

+  — ) 
m 


+  s  in  — ^-^'^)=m.  sinC 

^  m 

(30a) 


When  ip  is  not  a  multiple  of  m, 


sin  0  +  sin  (0  +  sin  (0  +  ••••  +  s  in  (0+11-^^^—^-^  r? )  =  q 

(30b) 

This  also  holds  true  for  a  cosine  sum  with  equivalent  arguments. 

Because  the  use  of  sines  and  cosines  in  K(t,r)  preserves 
symmetry  about  the  origin,  g  and  v  are  both  odd  integers. 
Therefore,  their  sums  and  differences  must  he  even. 
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With  ip  analogous  to  (g  +  v)  or  (g  -  v)  in  Equations  C30a) 
and  (30b) 


00  00 


MME(x,t)  =  ^  I  [t  a  sin(g6ot+^^)  +  j  A  sinCgcot-^)] 


g=l  V=1 


VTT.X  - 

P 
(31) 

where  for  the  first  sine  term  (g+v)  is  an  even  multiple  of  m 
and  for  the  second  sine  term  (g-v)  is  an  even  multiple  of  m. 

When  the  cosine  term  in  K(t,r)  is  not  zero,  the  traveling 
waves  may  also  be  found.  When 


B  =  -  qs  I  b  sin  (vtt/2)  K,  /v 
gv  TT  ^  max  g  V  /  y 


(32) 


where 


B  is  a  constant  whose  value  depends  on  g  and  v 

gv 


the  cosine  term  in  Equation  (21)  leads  to: 

00  00  m- 1 

MMF(x,t)=  I  I  I  B  cos  g(ajt  +  r27r/m)  cos  v  [  (ttx/T  )  +  (r2iT/m)  ] 
g=l  v=l  r=0  P 

(33) 

The  trigonometric  identity 

cos(a)  cos(S)  =  j  cos(a+B)  +  j  cos(a-3)  (34) 


yields : 


00  oo 


MMF(x,t)=  I  I  [j  B  cos(ga)t+  ^)  +  ^  B  cos  (ga:t  -  -^)] 
g= lv=l°  p  ^  p 


VTi^, 

P 
(35) 


where  for  the  first  cosine  term  (g  +v)  is  an  even  multiple  of  m 
and  for  the  second  cosine  term  (g  j  an  even  multiple  of  m. 
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Therefore,  the  complete  expression  for  the  time  variation 
of  the  magnetic  field  (armature  reaction  MMF  only)  referenced 
to  a  fixed  stator  due  to  the  specified  m^^^  internal  phase 
armature  current  is: 


MMF(x,t)=  I 

g  =  l 


I 

v  =  l 


It  a 


TT  X 

sin(ga)t  +  v— )  + 


y  B 


gV 


cos  (gujt  +  v— ) 


.  I''  ■  r  ^  '"X-.  m  n  r  ^  fX^T 

2  sin(gwt-v— )  -  j  ®av  <=05  (gwt  -  v— )  ] 

®  p  °  ‘  p 

(36) 

As  before,  for  the  first  two  terms,  (g  +v)  is  an  even  multiple 
of  m  and  for  the  last  two  terms,  (g  -v)  is  an  even  multiple 
of  m. 


Example 

See  Appendix  I  (Generator  Design)  for  a  reference  on 
the  generator  specifications  used  in  this  example. 

For  m=6,  the  terms  of  Equation  (36)  are  non-zero  when 
(g+v)  or  (g-v)  is  a  multiple  of  6.  Therefore, 

V  =  +  6k  -  g  (37) 


where 

k  is  an  integer  constant  equaling  0,  +1,  +2,  ... 
Table  6  gives  the  space  harmonic  indexes  (values  of  v) 
corresponding  to  the  odd  harmonics  of  the  time  variation  of 
the  armature  reaction  MMF.  A  forward  traveling  wave  is 


43 


represented  by  sin(gu)t  -  v— )  .  The  velocity  of  the  traveling 

P 

wave  components  relative  to  that  of  the  rotor  is  equal  to 
the  negative  time  harmonic  index  (-g)  divided  by  the  space 
harmonic  index  (v) . 


Table  6  Space  Harmonic  Indexes  for  Armature 
Reaction  Traveling  Waves 


The  traveling  wave  for  the  rotor  MMF  is  expressed  in 
Equation  (25),  The  harmonic  index  (v)  in  each  component  is 
the  same  for  the  time  and  space  variation. 

The  cancellation  of  two  orthogonal  functions  (Ref  10:188) 

A  sin(at)  +  B  sin(6t)  =  0  (38) 

requires  that  the  magnitudes  of  a  and  B  are  equal  or  that 
both  A  and  B  are  zero.  This  applies  to  the  cancellation  of 
armature  reaction  MMF  terms  (Equation  56)  and  rotor  MMF  terms 
Equation  (25). 


If  cancellation  is  to  apply  to  all  positions  (x)  and  all 
times  (t),  the  rotor  MMF  can  only  be  used  in  cancellation  of 
armature  reaction  MMF  terms  in  which  g  =  v.  In  otlier  words, 

rotor  MMF  will  only  cancel  armature  reaction  MMF  that  has  the 
same  relative  velocity  (g/v)  ;  the  velocity  of  the  rotor. 

Specification  of  rotor  MMF  is  not  sufficient  to  control 
the  generator  MMF  for  all  of  the  air  gap  at  all  times.  The 
armature  reaction  MMF  produces  both  space  and  time  variations 
that  are  not  cancelled.  This  limited  cancellation  is  more 
effectively  accomplished  by  the  specification  of  other  para¬ 
meters  such  as  pole  arc  (see  Section  II). 

Another  concern  with  this  approach  is  that  the  reluctance 
function  controls  the  flux  density  according  to  F.quation  (7). 
Thus,  even  if  a  certain  MMF  harmonic  is  cancelled,  there  must 
be  protection  from  the  reintroduction  of  that  harmonic  through 
the  variation  of  the  reluctance.  Therefore,  the  use  of  rotor 
MMF  is  not  a  favored  approach  in  the  control  of  generator 
harmonics . 
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VIII.  Reluctance  Harmonic  Control 

This  section  describes  an  attempt  to  control  phase  voltage 
harmonics  by  using  the  reluctance  to  control  tlie  flux  density 
harmonics . 

Neglecting  saturation  and  considering  the  air  gap  to  be 
dominant,  the  reluctance  along  the  air  gap  is  proportional  to 
the  air  gap  distance.  In  addition,  because  the  air  gap  rotates 
with  the  rotor,  the  reluctance  versus  distance  has  the  same 
shape  as  the  reluctance  versus  time  at  a  fixed  point  on  the 
stator.  Therefore,  although  the  reluctance  requirements  ma> 
be  solved  on  a  time  basis,  it  can  be  directly  related  back  to 
the  air  gap  and  thus  rotor  geometry.  MMF(x,t)  and  MMFft)  refer 
in  this  section  to  total  MMF. 

Flux  Density  from  MMF 

Just  as  it  is  the  combined  voltage  of  all  the  conductors 
connected  in  series  in  a  coil  group  which  is  specified  as 
phase  voltage,  it  is  the  sum  of  the  flux  densities  cutting 
these  same  conductors  that  must  be  controlled.  The  flux 
density  at  each  stator  slot  is  determined  by  the  MMl-  available 
to  force  the  flux  through  the  reluctance  at  that  point. 

The  MMF(x,t)  around  the  stator  is  known.  It  is  the 
superposition  of  the  MMF(x,t)  due  to  the  armature  reaction 
MMF  (Section  IV)  and  the  MMF  ^  due  to  the  rotor  coils 

(Section  V).  The  time  variation  MMF(t)  at  any  location  along 
the  stator  may  be  olitained  by  substituting  that  li'catian's 


46 


coordinates  into  the  MMl'(x,t)  expression.  For  a  j^iven  MMl  (t), 
the  value  of  reluctance  determines  the  magnetic  flux  density. 
I'herefore,  for  a  given  MMF(x,t),  the  desired  phase  voltage 
[Section  III)  may  be  liad  l)y  co:itrolling  tlie  reluctance. 

Tlie  reluctance  at  any  location  along  tlic  stator  periidiery 
is  also  a  function  of  time.  Therefore,  the  total  magnetic 
flux  dens ity  is: 


1  1  MMI'.It) 

'■TOT=  R^(t)  •  Area 
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K  is  a  proport i ona 1 i l y  conslant 
then  the  flux  density,  FApiation  (H),  is  a  constant  over  tliat 
pe  r i od  . 
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B  =  m/K 


(41) 


If  the  value  of  the  prof^ort  ional  ity  constant  K  is 
changed,  the  magnitude  of  the  flux  density  is  also  changed. 

The  basis  for  this  concept  involves  a  series  of  these  changes 
combined  with  careful  control  of  the  flux  density  harmonics. 

Reluctance  Specification  Concept 

The  magnetic  circuit  reluctance  function  is  periodic.  The 

reluctance  at  any  slot  is  the  same  function  as  the  reluctance 

at  a  previous  slot  except  for  a  time  delay.  If  the  MiMF(t) 

expression  for  each  slot  is  shifted  instead,  the  sum  on  the 

right  side  of  Equation  (39)  has  a  common  denominator.  Using 

* 

the  shifted  and  summed  MMF  (t)  reflected  into  the  first  slot, 
the  solution  of  reluctance  versus  time  will  be  simplified. 

A  pulse  width  current  modulation  technique  for  control 
of  line  voltage  harmonics  was  presented  by  Kr i shnamurthy  (Ref 
11:1269).  His  technique  is  reinterpreted  for  the  time  varia¬ 
tion  of  the  magnetic  circuit.  Here  permeance  is  modulated  to 
control  the  magnetic  field. 

The  resulting  flux  density  variation  B(t)  will  consist  of 
a  set  of  symmetrically  shaped  increments  for  each  set  of  har- 
monic(s)  controlled  (see  Figure  8c).  The  set  of  increment 
centers  for  the  n^^^  controlled  harmonic  is  given  the  name 
(see  Figure  8b)  which  corresponds  to  in  Reference  11. 
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(42) 


B  (C  +  At)  =  B  (C  -  At) 
V  n  V  ^  n  ' 


and 


B  (2k-rr  +  C  )  =  -B  (2kTT-C  )  =  B  [(2k+l)TT-C  ] 
V  n  V  n  n 


where 


-B^[  (2k  +  l)’i^C^] 
(43) 


k  is  an  integer  constant  equaling  0,  +1,  +2,  ... 
The  product  of  the  permeance  and  the  MMF  will  be  constant  in 
magnitude  over  the  increments,  but  will  change  sign  every  a 
radians . 

The  coefficient  of  the  n^^  harmonic  of  the  flux  density 
is 

8  B  ^  c 

n=l 


where 

d  is  the  width  of  an  increment  (maximum  At) 

c  is  the  number  of  harmonics  controlled  by  the 
increments  of  concern 

When  the  C  's  are  chosen  such  that: 
n 


c 

y  sin  vC  =0  (45) 

n=l 


by  can  be  made  equal  to  zero  for  c  different  harmonics.  There 
will  be  4c  pulses  per  2tt  electrical  radians. 
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Methodology 

The  harmonics  of  the  flux  density  are  controlled  by 
re  spec i f i cat  ion  of  the  air  gap  length  under  the  stator  slots 
during  certain  periods  .  The  T^  periods  correspond  to  the 
non- zero  periods  of  the  modulated  wave  claculated  for  no 
disturbance  of  the  n  previously  controlled  harmonics.  The 
air  gap  length  over  each  period  is  such  that  although  there 
is  variation  in  the  MMF,  the  flux  density  is  constant.  The 
process  is  accomplished  one  harmonic  at  a  time  with  the  funda¬ 
mental  completed  last. 

For  the  first  non- fundamental  harmonic  to  be  controlled, 
Uj ,  the  periods  T^  correspond  to  the  entire  pole  arc. 

Equation  (8)  is  rearranged  to  solve  for  the  length  of  the 
air  gap  over  the  desired  periods. 

u  MMF(t)„ 

6(t)  =  - - - (46) 


The  desired  values  for  the  flux  density  are  given  by  rearrange¬ 
ment  of  Equation  (10). 


’v  "  VI 


(47) 


For  the  second  harmonic  to  be  controlled,  V2,  the  air  gap 

is  respecified  over  the  periods  T. .  The  centers,  C  ,  of 

'■  1 

those  periods  are  chosen  in  each  fundamental  frequency  quarter 
cycle  so  that  Equation  (45)  is  satisfied  (c=l).  Therefore, 
there  will  be  no  contribution  to  tlie  Vj 


harmonic  from  tlie 


respecification.  The  width  of  the  periods  is  extended  to  the 
point  that  they  meet  each  other  or  the  edge  of  the  rotor  pole. 
The  respec i f i ca t ion  (addition  of  an  increment  of  permeance 
over  the  periods  T^]  is  done  so  that  the  time  harmonics 
of  the  flux  density  are  the  proper  value. 

Equation  (41)  is  used  to  calculate  the  flux  density  com¬ 
ponent  for  Vt  above  the  pole  arc  outside  of  T^ .  The  remainder 
of  the  necessary  flux  density,  A ,  must  come  from  the  period 
Tj .  The  flux  density  to  be  added  over  the  periods  T^^  must  be 
a  constant  according  to  the  modulation  technique. 

Reluctance  is  proportional  to  6 .  If  the  needed  reluctance 
is  6'  and  the  present  reluctance  is  6,  the  reluctance  to  be 
added  is  the  difference  Returning  to  the  distance  domain, 
the  new  air  gap  length  must  be  6'*  5/(6-6').  The  respecifi¬ 
cation  for  AB  becomes 


V 

If  T^  extends  over  the  boundaries  of  previous  T, Equation 
(48)  will  be  valid  as  long  as  the  proper  5  versus  time  is  used. 
Further  harmonics  are  controlled  by  additional  increments  of 
re luctance . 


cample 


See  Appendix  I  (Generator  Design)  for  a  reference  on  the 
generator  specifications  used  in  this  example. 
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The  third  harmonic  is  chosen  as  the  first  to  be  controlled. 
The  desired  flux  density  is  cero.  From  liquation  (46),  :  is 
infinite. 

is  now  determined  for  control  of  the  fifth  harmonic 
(.the  second  to  be  controlled)  without  disturbance  of  the  third 
harmonic.  The  sum 

1 

r 

^  ,  s  i  n  V  C 
n=  1  n  n 

is  zero  for  C  =  t:/3.  The  centers,  C,  ,  ri/3  radians  awav  from 
each  multiple  of  tt  radians  are  Tr/5,  2a/3,  4’t/5,  and  Sz/S. 

The  width  of  each  increment  may  be  up  to  tt/3  before  they  touch. 
At  this  full  width,  the  outside  limits  of  the  increments  T^^ 
are  tt/6  to  5tt/6  and  7tt/6  to  The  merged  increments  arc 

2t7/3  wide  corresponding  to  the  width  of  the  chosen  2z/3  pole 
arc.  In  other  words,  a  2tt/5  pole  arc  constant  air  gap  pro¬ 
duces  no  third  harmonic.  This  is  consistent  with  conventional 
pole  arc  design  (Ref  9:34). 

The  desired  flux  density  for  the  fifth  harmonic  comes  from 
Section  III  and  Fquation  (47).  All  of  this  flux  must  come  from. 
Tj^  because  the  rotor  pole  arc  docs  not  extend  past  T^^.  From 
Equation  (48)  and  the  fact  that  6  is  infinite,  6',  (the  limit 
of  [5 ' • 5/ (5 -6  '  ) ]  as  5  goes  to  infinity)  may  be  found. 

Next,  T^  is  determined  for  control  of  the  next  harmonic 
without  disturbance  of  the  third  and  fifth  harmonics.  ITpiation 
(45)  is  simultaneously  solved  for  v  =  3  and  v=5  yielding  Ci,=  .23,'i- 
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and  -  .433i'.  The  centci's  ol  tlie  B  (t)  pulses  at’e  .dsT':, 

^'2 

.433s^  .5()()a,  1.233",  1.433a,  l,5(){>n,  niul  1.7(j0': 

radians.  Tlie  periods  meoi  and  exterui  to  the  edite  oT  tlie  pole 
when  their  wiiltli  is  .HUvii  laidians.  Tlie  reluctance  and  I'lux 
density  waves  built  with  these  pei'iods  are  sliown  in  l  iyures 
8a  and  8c.  'I'he  a  i  i"  yaii  lenj^tli  over  these  periods  is  adjusted 
to  provide  the  necessary  flux  density  not  jirovided  outside  0  i' 
IT,.  I'he  flux  density  adjustment  continues  in  this  manner  for 
other  hai'/iionics.  'llie  reluctance  over  the  final  i)erjods  I 
which  affect  none  of  tlie  jireviously  controlled  harmonics  is 
adjustetl  for  the  fundamental  component. 

A  graphical  interpretation  lieljis  clarify  the  jirocess. 
i'lie  periods  of  non- infinite  reluctance,  shown  in  figure  8h, 
are  calculated  as  previously  described.  The  air  yap  lenyth 
(proportional  to  reluctance!  is  obtained  during  each  period 
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respectively.  The  curved  bound  of  the  shaded  areas  in  Figures 
8e  and  8f  represent  the  integrand  of  Equation  (50)  when  v=3  and 
5  respect i v c 1 y . 

The  signed  shaded  areas  in  Figures  8d  tlirough  8g  represent 
the  integrals  from  liquations  (49)  and  (50).  Algebraically,  the 
signed  areas  of  Figures  8d  through  8g  sum  to  zero.  This 
illustrates  that  Equations  (49)  and  (50)  are  zero  for  v=3  and 
for  v=5.  Thus,  the  harmonic  components  are  zero  for  these 
harmonic  '  because  of  the  chosen  reluctance  values  and  periods. 
The  illustration  is  similiar  for  other  higher  harmonics. 

How  Man >'  Harmonics 

The  process  as  described,  approaches  the  desired  flu.x 
density  wave  shape  one  harmonic  at  a  time.  To  see  the 
advantage  of  this  approach,  consider  specification  of  the 
reluctance  as  a  composite  vvave  shape. 

5(t)  = 

where 


5(t)  is  a  time  varying  air  gap 

This  requires  a  specific  finite  reluctance  at  all  times.  When 
the  large  air  gap  betvveen  the  poles  is  under  the  coil  group, 
the  reluctance  is  not  variable.  This  shortcoming  affects  all 
the  harmonics.  Therefore,  the  chosen  method,  because  the 
important  harmonics  can  be  directly  controlled  has  an  advantage. 
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The  number  of  harmonics  uTiich  can  be  controlled  is 
dependent  on  those  chosen  and  tlie  pole  arc.  Equation  (45) 
yields  =  .185tt,  =  .3447t,  and  C,  =  .  4337t,  for  control 

of  the  3rd,  5th,  and  7th  harmonics.  If  the  pole  arc  length 
was  less  than  .633t:,  C,  would  be  at  the  pole  edge  and  control 
of  these  three  harmonics  could  not  be  achieved.  Therefore, 
the  harmonics  and  the  pole  arc  must  be  selected  consistent 
with  the  methodology. 

The  Contradiction 

The  concept  of  harmonic  control  appears  to  be  practical 

once  MMF„  (the  reflected  MMF)  is  obtained  and  the  control 
K 

harmonics  are  selected.  However,  obtaining  MMFj^  is  not  an 
easy  problem. 

The  proposed  shift  of  MMF  at  consecutive  slots  1  and  2 
is 

MMF^*(t)  =  MMF2(t-At)  (52) 

where 

MMFj^  (t)  is  the  shifted  MMF 

The  reluctance  at  one  slot  is  the  same  function  as  the 
reluctance  at  the  previous  slot  except  for  a  delay  of  At. 

In  other  words 

Rj(t)  =  R2(t  '  At)  (53) 

An  expression  for  the  flux  density  at  the  second  slot  is 
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(t  -  At) 


(54) 


MMF^  (t-At) 


The  shifted  MMF  and  reluctance  arc  then  plugged  into 
Equation  (54).  The  contradiction  is  that  the  flux  density 
B2(t-At)  reflected  into  the  first  slot  is  not  synchronized 
with  the  flux  density  at  the  first  slot  Bj^(t).  Such  a  sum 
makes  little  sense. 

A  general  solution  for  a  reflected  B2(t)  synchronized 
with  Bj^(t)  and  with  the  knowledge  of  MMF2(t),  but  not  R2(t) 
exceeds  the  resources  of  this  thesis  effort.  Therefore, 
another  approach  is  preferred. 
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IX. 


I  nc  rciiicnt  a  I  Pcluctaricc  Spec  i  f  i  c a  t  i  on 


} 

This  section  (.Icscrihes  th.c  final  iiictbotl  in  uhicli  the 
rotor  geometry  specification  is  analyzed  for  control  of  tlie 
generated  voltage  wave  shajic.  The  reluctance  o  1'  the  air 
gap  is  specified  as  tfte  traveling  rotor  is  positioned  under 
consecutive  stator  slots.  .lust  as  in  Section  I'lII,  reluc¬ 
tance  versus  time  may  also  he  related  to  air  gap  distance 
t’ersus  time  and  tluis  to  rotor  geometr\-.  The  s  jtec  i  f  i  c  a  t  i  on 
j  jn'ogresses  with,  time  along  with  the  rctiuired  wave  shape. 

Tlierefore,  individual  harmonic  components  arc  not  of  concern. 

to  I' lux  Itensity 

The  total  versus  time  ivliich  was  developed  in  .Section 

I’ll  is  a  useful  starting  jilace  for  this  develojniUMit  .  The 
reluctance  must  be  specified  so  that  the  resulting  magnetic 
flux  density  will  t'e  in  accordance  with  the  magnetic  flux 
density  required,  as  was  discussed  in  .Section  111.  .X  i  r  gap 
dominance,  the  lack  of  fringing,  and  tlie  finite  i^ermeance 
effects  arc  assumed  lierc  as  they  were  in  Section  1  !  I  . 

The  reluctance  at  any  slots  is  th.e  same  function  as  the 
reluctance  at  the  previous  slot  except  for  a  time  delay  of  .dt. 

(  > 
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Therefore,  the  total  flux  density  may  be  expressed  as: 


where 


B^QjCt) 
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B  Ct) 
n  ^ 


^  MMF  ( t ) 

L  - 

n=l  Rj (t -nAt ) 


n  is  the  individual  stator  slot  number 
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Reluctance  Specification 

When  the  gap  between  salient  poles  is  under  a  slot,  the 
reluctance  at  that  slot  may  be  approximated  as  infinite.  When 
the  rotor  turns  so  that  its  edge  is  between  the  first  and 
second  conductor  of  the  relevant  coil  group.  Equation  (S6) 
has  a  non-zero  term  only  for  n=l.  The  desired  increment  of 
flux  density,  AB,  and  the  available  MME  determine  the  reluc¬ 
tance  required  while  the  rotor  edge  is  between  slots  1  and  2. 

MMF  (t) 

After  a  delay  of  t,  that  same  reluctance  will  appear 
under  the  second  slot.  The  MMF^  at  (t+At)  and  the  already 
specified  reluctance  will  determine  the  flux  density  at  the 
second  slot.  The  remainder  of  the  desired  total  flux  density 
must  come  under  the  first  slot.  By  using  MMF(t+At)  and  the 
residual  requirement  of  reluctance,  AB,  in  Equation  (B7),  the 
reluctance  at  the  next  increment  may  be  determined.  The 
reluctance  specification  process  is  repeated  as  the  rotor 
progesses  under  the  coil  group  and  is  completed  once  the  new 
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rotor  position  is  tt  radians  from  whore  it  was  when  the 
specification  began. 

Ihe  specification  of  tlie  reluctance  using  t  ii  i  s  method  is 
very  simple.  However,  insight  is  required  in  the  choice  of 
the  phase  voltage  wave  sha[ie.  After  tlie  trailing  edge  of  the 
rotor  pole  face  passes  the  first  slot,  tiie  reluctance  specifi¬ 
cation  is  complete.  The  phase  voltage  and  the  recjuired  mag¬ 
netic  flu.x  densit)'  must  come  from  the  remaining  MMf  tiiiie 
variation  and  the  reluctance  specified  to  this  point. 

lixampl  e 

See  Appendix  I  ffenerator  Design)  for  a  reference  on  tlie 
generator  spec  i  f i ca t i ons  used  in  this  example. 

The  proof  that  not  just  any  phase  voltage  wave  sha[K'  will 
work  is  seen  by  examination  of  the  example  developed  so  far. 

Table  5  gives  the  time  variation  of  the  total  MMT  at  each  slot 
for  the  Figure  5  voltage  variation.  At  a/t  =  -5r/l5,  none  of 
the  pole  arc  is  under  the  phase  C  coil  group.  At  '-t  =  -  4'i  / 1  5  , 
the  pole  will  be  under  the  first  conductor  of  the  coil  group. 

The  MMF  at  v.t  =  -4Ti/15  is  negative.  Additional  rotor  MMF 
could  be  added  here  in  order  to  ensure  a  positive  magnetic 
flux  density  with  a  realicable  reluctance.  However,  there  is 
another  difficulty.  Tlie  desired  magnetic  flux  density  and 
therefore  the  desired  voltage  wave  shape  is  flat  between  .  t  -  2  •■  /  1  .S 
and  wt  =  +  /  1 -S  .  At  <.ut  =  -  2tt/ 1  .S  ,  a  certain  combination  of  reluc¬ 

tances  will  be  under  slots  1  ,  2,  and  a  and  w  i  t  ti  t  lu-  < 
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(Table  5  entries)  will  produce  the  desired  magnetic  flux 
density.  However,  at  (i)t=-7i/15,  the  reluctances  are  exposed 
to  greater  MMFs.  The  new  increment  under  slot  1  also  sees 
a  positive  MMF .  Then  the  reluctances  wliich  provide  the  correct 
amount  of  flux  density  at  a)t=-27r/15  arc  bound  to  provide  too 
much  at  ojt=-7T/15.  Therefore,  this  choice  of  a  phase  voltage 
wave  shape  is  not  consistent  for  wave  shape  control  using 
only  rotor  geometry. 

Triangular  Phase  Voltage  Wave  Shape  Example 

The  flexibility  of  the  phase  voltage  wave  shape  for  a 
required  load  voltage  wave  shape  as  discussed  in  Section  III 
is  valuable.  Figure  9  shows  a  phase  voltage  wave  shape  wliicli 
produces  a  DC  load  voltage.  Notice  that  it  has  an  upward 
ramping  voltage  as  the  number  of  slots  with  finite  reluctances 
increases  and  a  downward  ramping  voltage  as  the  number  of  slots 
with  finite  reluctances  decreases.  In  this  example,  the  peak 
phase  voltage  equals  the  load  voltage  and  the  pole  arc  is  only 
tt/3  electrical  radians.  Tables  7,  8,  and  9  display  the  voltage, 
the  magnetic  flux  density,  the  current,  and  the  MMF  variation 
at  each  slot  versus  time. 

The  reluctance  for  this  new  wave  shape  is  introduced 
under  the  first  conductor  at  'jjt=  -  So/SO .  The  rotor  pole  is 
centered  under  phase  C  at  ojt  =  -5TT/50.  It  is  under  only  one 
conductor  at  jt  =  137T/30  as  it  is  leaving  the  coil  group. 
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Table  10  shows  the  progressive  calculation  of  the  reluc¬ 
tances  required  for  the  Figure  9a  triangular  voltage  wave  shape. 
Notice  that  Table  10  stops  with  the  determination  that  a  nega¬ 
tive  reluctance  is  required.  A  rotor  MMF  of  1400  ampere-turns 
was  chosen  instead  of  the  533  ampere-turns  used  prev'iously  in 
order  to  prevent  the  calculation  of  a  negative  reluctance  for 
the  very  first  entry.  This  wave  shape  choice  is  not  appropriate 
for  rotor  geometry  control. 

From  the  shaded  area  of  Table  9,  it  is  observed  that  an 
equation  for  the  magnetic  flux  density  (similar  to  r:quation 
39)  may  be  written  for  each  line  (at  each  time  increment). 

Nine  equations  with  the  specified  MMF  determine  the  required 
reluctances.  Yet,  only  five  reluctances  are  variable.  The 
armature  current  and  rotor  MMF  provide  two  additional  variables. 
It  is  clear  that  design  using  the  method  of  this  section  is  a 
form  of  over- spec i f icat ion .  There  might  be  a  waveform  that 
fits  all  of  the  equations  and  provides  the  desired  load  vol¬ 
tage,  but  none  were  found. 


Table  7  Normalized  Phase  Current  and  Voltage 
(Triangle  Wave  Shape)  (Figure  9a) 


TIME 

(wt) 

VOLTAGE 

CURRENT 

A 

B 

C 

A' 

B' 

C 

-  Stt/jO 

0 

1 

0 

-1 

0 

1 
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0 

.8 

-.2 

-1 

.36 

.  96 
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0 

.6 

-  .4 

-1 

.  64 

.  84 
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0 

.  4 

-.6 

-1 

.  84 

.64 
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.2 
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-1 
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57r/30 

0 

0 

1 

- 1 

1 

0 
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7 

0 
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1 
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.4 

0 
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-  .84 

1 

-  .64 
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.6 

0 

-  .4 
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1 
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.  8 

0 

-  ^  2 

-  .34 

1 

-  .96 
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1 

0 

0 

0 

1 

-1 
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.8 
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0 

.  34 

.  96 

-1 
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.6 
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0 
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-1 
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.4 
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0 

.  84 
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,  2 

-  .8 

0 
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-1 
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1 

0 
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Table  8  Normalized  Armature  Reaction  MMF  Phase 
C  Slots  (Triangle  Wav^e  Shape) 


TIMe\  X 
(‘"t)  \ 

77/30 

377/30 

577/30 

7  77/3  0 

977/3  0 

-Stt/SO 

-1.80 

-1.40 
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-  .  60 

-  .  20 

-37r/30 

-2.13 

-1.74 

-1.36 

-  .  98 

-  .  59 

-Tr/30 
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-1 . 98 
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71/30 
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-1.67 
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1 

00 
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1177/30 

-.97 
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-1.98 

-2.31 

1  377/30 

-.57 

-  .96 

-1.34 

-1.72 

-2.11 

1  577/30 

-  .  20 

-  .60 

-1.00 

-1.40 

-1.80 

1777/30 

.  18 

-.22 

-  .  62 

-1.02 

-1.42 

1  977/  30 

.60 

.20 

-.20 

-  .  60 

-1.00 

2177/30 

1 . 00 

.60 

.  20 

-  .  20 

-  .  60 

2  377/30 

1.40 

1.00 

.60 

2i 

-  .  20 

2577/30 

1.80 

1.40 

1 . 00 

.  oO 

.  20 
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Table  9  Total  MMF  at  Phase  C  Slots  (Triangle  Wave  Shape! 
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Figure  9  Triangle  I’liase  Voltage  ani.1  (Tirrent  Wave 
S  li  a  p  e  s 


X .  Other  Magnetic  nfrpct? 


This  section  discusses  the  effects  of  flux  leakage, 
flux  fringing,  and  flux  saturation  on  the  desired  generated 
\-oltagc  wave  shape.  Finite  Element  analysis,  an  important 
analytical  tool,  is  introduced  as  a  computer-aided  design 
method  used  in  fine  tuning  and  checking  the  generator 
design.  One  example  of  a  magnetic  finite  element  analysis 
program  is  introduced. 

Introduction 

.A,ll  of  the  procedures,  methods,  and  examples  in  this 
thesis  effort  assume  zero  flux  leakage,  no  flux  fringing, 
and  no  saturation  in  the  generator's  magnetic  structure. 
These  are  important  assumptions  because  calculating  flux 
leakage,  flux  fringing,  and  saturation  require  a  knowledge 
of  the  complete  magnetic  structure;  something  that  is  not 
available  in  the  initial  design  stages.  Generator  design 
is  by  nature  an  iterative  process  due  to  these  three 
problems . 

When  a  generator  is  designed  for  a  specific  load,  the 
system  inductance  is  an  important  specification.  The  sys 
tern  inductance  depends  on  the  load  inductance,  the  trans¬ 
mission  lino  inductance,  and  also  tlie  generator  inductance. 
The  generator  inductance  can  easily  be  the  most  significant 
of  the  three  inductances.  The  generator  inductance  depends 


on  the  geometry  of  the  magnetic  circuit.  Designing  the 
magnetic  circuit,  however,  depends  on  a  know'lcdge  of  the 


system  inductance.  Normally,  a  "first  cut"  design  is  made 
from  which  a  "first  cut"  system  inductance  can  he  calcu¬ 
lated.  Using  th.is  inductance,  a  "second  cut"  design  is 
made  from  which  a  "second  cut"  system  inductance  can  be 
calculated.  This  process  continues  until  it  converges  on 
a  final  design. 

This  thesis  effort  is  concerned  only  with  the  "first 
cut"  design.  The  rest  of  the  design  procedure  is  left  for 
further  study. 

Leakage  Flux 

In  the  generator,  the  rotor  magnetic  flux  is  generated 
in  the  rotor  coils.  The  flux  travels  through  the  rotor 
poles  and  into  the  stator  magne+ic  structure  cutting  the 
armature  coils  and  th.us,  generating  the  output  voltage. 

In  a  salient  pole  generator,  there  is  a  large  air  gap 
between  two  adjacent  poles.  The  magnetic  reluctance  is 
very  high,  but  '  have  a  finite  value.  Since  there  is 
an  MMF  drop  betwt  two  adjacent  poles,  there  will  be  mag¬ 
netic  flux  crossing  the  gap  between  the  poles  bypassing 
the  stator.  All  of  thl?  leakage  flux  which  bypasses  the 
stator  does  not  l.ave  a  significant  impact  on  the  shape  of 
the  desired  outfut  voltage  waveform,  although  it  docs 
increase  the  gene!'it<'r  inductance,  X  .  To  compensate  for 


this  leakage  flux,  all  the  designer  needs  to  do  is  increase 
the  numher  of  ampere  turns  in  the  rotor  coils. 

Fringing  Flux 

Another  important  assumption  made  in  this  thesis 
effort  is  that  all  the  flux  generated  in  the  rotor  passes 
through  the  air  gap,  past  the  armature  coils  and  into  the 
stator  core  in  a  totally  radial  direction.  This  is  an 
important  "first  cut”  assumption  and  can  be  corrected  for 
in  further  design  iterations. 

In  a  realistic  generator,  the  majority  of  the  magnetic 
flux  does  go  in  the  radial  direction.  However,  there  is 
a  significant  amount  of  flux  which  does  not.  This  flux 
crosses  the  air  gap  in  directions  defined  by  the  magnetic 
flux  density  vector  ?  causing  a  distortion  of  the  desired 
output  voltage  wave  shape.  The  flux  can  also  short  cut 
across  several  armature  slots  instead  of  going  straight 
through . 

.Accounting  and  correcting  for  this  fringing  flux  so 
that  the  correct  flux  links  the  armature  coils  is  no  simple 
matter.  Where  and  how  the  flux  fringes  depends  on  the 
shape  of  the  generator's  magnetic  structure.  The  correc¬ 
tion  for  the  flux  fringing  is  therefore  an  iterative  design 


process . 


Magnetic  Saturation 

The  last  important  assumption  made  in  this  thesis 
effort  is  that  the  relationship  between  the  MMF,  the 
reluctance  of  the  magnetic  material,  and  the  resultant 
flux  density  in  the  magnetic  material  is  linear.  This 
last  assumption  said  another  way  means  that  the  analysis 
in  this  thesis  effort  ignores  the  effects  of  saturation  of 
the  magnetic  material. 

In  a  realistic  generator,  this  relationship  is  defi¬ 
nitely  non-linear.  When  the  ^fMF  dropped  across  the  magne¬ 
tic  structure  is  a  given  value,  a  certain  amount  of 
magnetic  flux  will  "flow"  through  the  material.  As  the 
NIMF  increases,  the  flux  flowing  through  the  material  also 
tends  to  increase.  However,  at  high  enough  values  of  flux 
density,  the  magnetic  material  saturates  and  more  and  more 
MMF  is  required  for  smaller  increases  in  the  flux. 

As  was  true  for  the  flux  fringing,  the  locations  and 
the  severity  of  the  flux  saturation  are  determined  by  the 
structure  of  the  magnetic  circuit.  Compensating  for  this 
problem  can  also  be  an  iterative  design  procedure. 

Finite  Element  Analysis 

The  calculations  involved  in  designing  a  generator 
can  be  exceedingly  complex  and  cumbersome.  This  fact  makes 
the  digital  computer  a  welcome  design  aid.  One  particular 
procedure  ideally  suited  for  computer-aided  design  is 
Finite  Flcment  analysis. 
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In  Finite  Element  analysis,  a  cross  section  of  the 
generator  is  taken  perpendicular  to  the  axle.  Then  the 
assumption  is  made  that  any  leakage  flux  or  fringing  flux 
is  in  the  plane  of  the  cross  section.  That  is,  there  is 
no  leakage  flux  or  fringing  flux  between  any  two  "adjacent, 
infinitely  thin"  cross  sections.  The  two  dimensional  cross 
section  is  then  divided  up  into  triangular  finite  elements. 
Each  element  is  restricted  to  one  particular  type  of  mate¬ 
rial  (steel,  air,  insulation,  ...).  During  the  computa¬ 
tions,  each  triangular  finite  element  is  considered  to 
have  a  constant  magnetic  flux  density  over  the  entire 
element.  See  (Ref  8)  and  (Ref  13)  for  further  information 
on  the  Finite  Element  method. 

MAGNETIC  is  a  computer  program  that  uses  the  Finite 
Element  analysis  technique  (Ref  4)  to  analyze  magnetic 
structures.  MAGNETIC  computes  the  magnetic  fields,  induced 
currents  and  losses,  and  the  forces  of  a  magnetic  device 
that  has  been  represented  as  a  collection  of  finite  ele¬ 
ments  . 


"MAGNETIC  aids  the  design  of  magnetic 
devices  because  accurate  prediction  of  devices 
performance  usually  requires  knowledge  of  the 
magnetic  field.  MAGNETIC  is  a  means  of  calcu¬ 
lating  the  field  that  allows  for  fringing, 
saturation,  permanent  magnets,  and  complicated 
shapes.  Because  the  program  includes  non-linear 
B-H  curves  and  computes  inductances,  forces,  and 
eddy  current  losses,  it  can  help  the  designer  to 
minimize  the  materials  used  in  a  device  and  maxi¬ 
mize  its  efficiency."  (Ref  4:11) 


Two  problems  in  generator  design,  flux  fringing  and 
saturation,  force  the  design  procedure  to  be  iterative. 
MAGNETIC  can  be  a  welcome  aid  in  correcting  for  flux 
fringing  and  saturation  as  well  as  flux  leakage.  After 
the  designer  has  completed  the  "first  cut"  design,  the 
design  can  he  analyzed  by  MAGNETIC.  MAGNETIC  shows  through 
its  output,  where  and  to  what  extent  flux  fringing,  satura¬ 
tion,  and  flux  leakage  are  occurring.  With  this  informa¬ 
tion,  the  designer  is  able  to  take  corrective  action  and 
come  up  with  a  "second  cut"  design.  This  procedure  con¬ 
tinues  until  an  acceptable  design  is  obtained. 

Tn  addition  to  tabular  listings  of  the  various  output 
parameters,  MAGNETIC  produces  a  Mesh  plot  and  a  Flux  plot. 
The  Mesh  plot  portrays  the  selected  generator  cross  section 
showing  all  the  triangular  finite  elements.  The  Flux  plot 
shows  the  lines  of  magnetic  flux  in  the  generator  cross 
section.  Figure  10  is  the  Mesh  plot  for  the  sample  gener¬ 


ator  Cl^eference  Appendix  I)  .  The  Flux  plot  of  the  sample 
generator  is  not  shown  because  that  portion  of  MAGNETIC 
currently  does  not  work. 


XI.  Conclusions  and  Recommendations 


Conclusions 

The  purpose  of  this  thesis  effort  was  to  investigate 
the  use  of  rotor  geometry  as  a  specification  in  the  control 
of  generator  voltage  and  flux  density  harmonics.  It  was 
found  that  armature  reaction  IIMF  could  be  expressed  as  an 
infinite  set  of  traveling  waves  with  an  infinite  set  of 
velocities  (Section  IV)  and  (Equation  23).  Because  the 
rotor  has  only  the  velocity  of  the  rotor,  it  cannot 

cancel  armature  reaction  everywhere  in  the  generator 
throughout  all  time  (Section  VII). 

Rotor  geometry  was  put  to  use  in  a  limited  control  of 
flux  density  harmonics  at  the  stator  slots.  The  reflection 
of  the  MMFs  present  at  all  of  the  slots  of  a  coil  group 
into  a  single  slot  allows  the  possible  calculation  of  an 
appropriate  rotor  geometry.  This  reflection  turned  out  to 
exceed  the  resources  of  this  thesis  effort  and  was  not 
accomplished  (Section  VIII). 

Specifying  the  air  gap  reluctance  incrementally  as 
the  rotor  moves  under  successive  stator  slots  was  also 
tried  as  a  means  of  harmonic  control.  This  method  was 
unworkable  for  all  wave  shapes  chosen  because  it  produced 
negative  NiMF  requirements  and  conflicting  geometry  speci- 
f icat ions . 
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Recommendations 


The  control  of  flux  density  is  an  important  consid¬ 
eration  when  working  to  improve  generator  efficiency  and 
performance.  Because  of  the  potential  that  flux  density 
control  has,  several  recommendations  for  further  study 
are  presented. 

Before  any  hardware  is  fabricated,  it  is  recommended 
that  the  designer  calculate  the  benefits  due  to  the  flux 
density  control  design.  Knowing  that  energy  is  conserved 
betw'een  the  input  and  the  output  of  the  generator,  the 
designer  will  be  able  to  calculate  the  design's  value. 

A  recommended  approach  for  reduction  of  losses  is  to 
investigate  flux  density  control  at  the  damper  bars  instead 
of  the  stator  windings.  This  investigation  however,  will 
not  control  the  generator  output  harmonics. 

A  different  approach  might  be  to  use  stator  geometry 
as  a  design  specification  in  harmonic  control.  This  will 
not  control  the  MMF  at  all  time  and  all  space  by  itself, 
but  may  provide  easier  control  of  the  harmonics  at  a  given 
location . 

The  method  in  Section  VII  could  be  further  investi¬ 
gated  if  a  method  were  derived  for  the  division  of  one 
infinite  series  with  progressive  harmonic  arguments  by 
another  such  scries.  The  Section  VIII  approacli  could  be 
further  investigated  in  terms  of  a  general  wave  shape 


(sec  the  headinj^  "The  (’out  rad  i  ct  i  on"  in  Section  VI  !1). 
liither  method  might  provide  a  workable  wave  sliapc  or  prove 
tl'.at  the  rotor  geometry  concept  is  in  general  imjiossihle. 

The  general  requirement  for  phase  voltage  wave  shapes 
tl'.at  produce  a  IK,'  output  Crom  a  tliree  ]ihase  full  wave 
!•'  r  i  d  g  e  rectifier  is: 
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The  final  recommendation  is  to  include  the  flux 
fringing,  flux  leakage,  and  magnetic  saturation  problems 
in  the  calculation  (only  if  the  calculations  are  successful 
without  these  three  problems)  for  the  desired  phase  voltage 
wave  shape.  It  should  be  recognized  that  the  problem  will 
now  he  iterative  and  Finite  Element  analysis  (Section  X) 
will  be  helpful. 
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Appendix  I  Sample  Generator 


Specifications  and  Dimensions  for 
the  Sample  Generator 
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The  development  of  the  rotor  geometry  modification 
example  in  this  thesis  requires  a  constant  air  gap,  salient 
pole,  synchronous  generator  to  begin  the  investigation. 

The  generator  selected  is  basically  the  one  designed  for 
tlie  final  exam  in  Air  Force  Institute  of  Technology  course 
EF  6.07  (Electromagnetic  Components  Design,  Spring,  1980). 

A  few  modifications  arc  made  due  to  thesis  requirements. 

The  design  procedure  for  the  sample  generator  is  not  pre¬ 
sented.  Specifications  for  the  sample  generator  are 
presented  in  Table  11.  A  cross  sectional  drawing  (per¬ 
pendicular  to  the  axle)  of  the  generator  is  presented  in 
Figure  11. 

Th ■  design  modifications  developed  in  this  thesis 
apply  equally  as  well  to  any  salient  pole,  synchronous 
generator  as  they  do  in  the  generator  used  to  begin  the 
invest igat ion . 

Notice  that  the  generator  design  is  in  inches.  -Ml 
units  shown  in  this  appendix  are  converted  to  the  MKS 
system  for  use  in  the  body  of  this  thesis. 
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Table  11 


Sample  Generator  Dimensions  and  Specifications 


Machine  Specifications 

WYE  connected  three  phase  output 

460  volts  WfS  line  to  neutral  phase  voltage 

72.464  amperes  RMS  phase  current 

Load  power:  100  kilowatts 

Frequency;  400  hertz 

No  rotor  or  stator  vents 

Straight  slot  stator  design 

60  slots 

4  poles 

0.015  inch  constant  air  gap  over  the  poles 

Synchronous  speed:  12000  revolutions  per  minute 

Rotor  tip  speed;  393  feet  per  second 

Rotor  MMF :  533  ampere  turns  per  pole 

Maximum  flux  density:  1  Tesla 

Pole  arc;  60  mechanical  degrees,  a  =  0.667 

Rotor  length;  25.347  inches  ^ 

Trancor  "T"  steel  with  0.007  inch  laminations 
Synchronous  salient  pole  design 


Winding  Specifications 

Six  internally  generated  phases,  m  =  6 
Armature  conductors  have  no  skew,  K  =  1.0 
q  =  5 

One  series  current  path  per  output  phase 

Full  pitch  stator  winding,.  W/T  =  1.0,  K  =  +  1.0 

Two  layer  lap  design  ^  P  • 

0.107  by  0.107  inch  square  conductors 

15000  amperes  per  square  inch  maximum  current  density 

80  volts  per  mil  minimum  insulation  thickness 

T  =15  slots 
P 


Assumptions  (for  "first  cut"  calculations) 


Zero  flux  leakage 
Zero  leakage  reactance 
PF  =  1.0 
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Figure  11  Cross  Section  of  the  Sample  Generator 
(Specifications  and  Dimensions) 
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Appendix  I 


RGRID  and  TRIA  Card  Generation 


Computer  Program  to  Generate  RGRID  and  TRIA  Cards 
Required  by  MAGNETIC  (Ref  4)  for  the  Air  Gap  and  Pol 
Shoe  Regions  of  the  Sample  Generator  (Appendix  I) 


Appendix  II  RGRIP  and  TRIA  Card  Heneration 


The  following  program  listing  generates  the  RGRID  and 
TRIA  cards  required  by  MAGNETIC  (Ref  4)  for  the  air  gap 
and  pole  shoe  regions  of  the  sample  generator  (Appendix  I). 

Refer  to  Figure  12  and  to  the  program  listing  for  an 
explanation  of  how  to  use  the  program.  Basically,  to 
modify  the  program  for  another  generator,  the  following 
must  be  accomplished: 

1)  Fill  in  the  appropriate  numbers  in  the  three  data 
statements . 

2)  Calculate  values  for  the  theta  and  radial  coordi¬ 
nates  for  the  points  A,  B,  C,  and  D. 

3)  Write  a  statement  function  that  describes  the 
surface  of  the  pole  shoe. 

The  internal  program  structure  is  not  discussed  because 
it  is  not  important  for  the  program’s  use. 
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Figure  12  Cross  Section  of  the  Sample  Generator  (RGRID  and  TRIA  Cards) 
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PROGRAM  THESISl  ( INPUT» OUTPUT , PUMGH,  TAPE£  =  0UTPUT ,T APP7=PU <CH) 

T“TS  PROGRAM  GEMEFATPS  ONLY  THE  AlR  3A»  AND  =0LE  SHOE  RGRTO 
AK'O  CAFDS  F.EOUI^ED  BY  THE  PROGRAM  MAGNETIC  (REE  <* )  FOR 

Tur  SAMOLE  GtNEFATCR 

V.flRIAOLE  OIMEUSIONS 

TNG  -  ANGLE  FOR  P^OIVIOUAL  RGRIO  CARDS 

tT£Y  -  KEYrI  DESIGNATES  A  TRIA  IS  AIR,  <£Y  =  '’  DESIGNATES  ‘^TEEL 
LNG  -  LNG+1  =  FIRST  GRID  NUMBER  ASSIGNED 

LNS  -  LNS+1  =  LOWEST  GuIO  NUM^Er;  ON  THE  STATDE  AIR  GAP  EOUNOARY 

LNT  -  LNT+l  r  FI'ST  TkIA  NUMBER  ASSIGNED 

IPG  '  LFG4-1  =  FIt-ST  GRIO  NUMBER  ON  THE  RB  ARS 

rr  -  mumper  of  c-GRIO  oCINTS  DN  segment  KF  (INCLUDE?  THE  ENDS) 

RPA  -  NUMBER  OF  PGRID  POINTS  ON  SEGMENT  KFA  (INCLUDES  THE  ENDS) 

Mcg  _  N'lMBEv  OF  RGRIO  POINTS  ON  SEGMENT  KFB  (INCLUDES  THE  ENDS) 

NG  -  RGRIO  NUNBER  TO  BE  PUNCHED 

NGl  -  NUMBER  OF  FGRID  POINTS  ON  SEGMENT  KGl  (INCLUDES  THE  ENDS) 

)‘P  -  NUMBFE  CF  FGRTD  POINTS  ON  SEGMENT  K®  (P'CLUnEE  THF  ENDS) 

^PA  -  NUMBE-  OF  RGRID  POINTS  ON  SEGMENT  KRA  (INCLUDES  THE  ENDS) 

‘■PB  -  NUMBEF  OF  RGRID  POINTS  ON  SEGMENT  KP^  (INCLUDE?  THE  ENDS) 

NT  ..  tbIA  HUMBER  TO  BE  PUNCHED 
-  DOMSTANT  3. 1<+Ii3926535 

RAO  -  RADIAL  CTSTANCE  TO  the  POLE  FACE  SURFACE 

BGl  -  RADIAL  distance  TO  THE  mioOLE  OF  THE  AIR  GAP 

YAB,  XCC,  XCD,  XT,  XY,  X7,  ANB  YCC  ARE  USED  TO  CALCULATE 
the  COO^<D:NATES  FCR  THE  POINTS  A,  3,  C,  AND  D. 

PA,  RD,  RC,  AND  PC  ARE  THE  RAPin  DODROINaTES  OF  A,  P,  C,  ANO  0. 
TA,  TB,  TC,  ANO  TP  ARE  THE  THETA  COORDINATES  OF  A,  B,  C,  AND  D, 

THE  following  variables  ARE  USEO  INTERNALLY  TN  THE  OROGRAM  ONLYt 
DKTB  ONFP  ONFPA  DNFPB  DNGFM?  DNSGl 

I  KP  KT  NFAMl  NFBM.l  NFMl 

NGIMI  NPAMl  NP9MI  N»M1 

the  FOLLOWING  ARE  DESIGNATED  LINE  SEGM'^NTS  -  SEE  FIGURE  12 

<F  KFA  KFB  <Gl  KP  <PA 

<P9 

SUBROUTINES  REOUIREOJ 

TRIAl  -  PUMCHE?  TRIA  CAROS  FOR  TY’E  I  FINITE  ELEMENTS 
TRIA2  -  PUNCHES  TRIA  CAROS  for  TYPE  2  FINITE  ELEMENTS 
SEE  triaI  and  TKIP2  FOR  FURTHER  EXPLANATION, 

DOUBLE  PRECISION  TA, TB, TC ,TD , RA , R3 , RD , PO 

OOUBLF  FRECISIDN  ANG,RAO,PI 

DOUBLE  PRECISICN  XY, X Z, XT  ,XCC , YCC, XCD , XA 3 
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1.>00^>C}  «  \.>  <•>  t>  k.>  t>  <■>(.>  t» 


OMft  NFr.,NF,NF'>,NPA,NP,NP9,»4Gl,NS  /II ,  1 , 1 1,  f,  2  3  ,61, 61/ 

TATA  LNG,LNT,LPG,LN9  /27  5  ,<*r.  t> ,  24 , 2 IV  / 

DATA  °G1  /3.;!.''5/ 


3  T‘'TS  '^TATF'^fkT  F'J».CT:0N  represents  the  SJ^FA;:  PF  i-qr  poLE 
:  IT  IS  inf  distance  versus  X  (HE  ORIGIN  IS  IN  TH'^ 

Q  C-f.‘T-5  OF  THE  PCLt)  A90VE  THE  PE  ARC.  IN  'THIS  CASE,  t  CONSTANT 
3  AT:='  gap  has  PEEK  SFECIFIEO.  THE  OnLT  RESTPICTION  IS  THA’  THE 
C  °CLE  PACE  nay  NOT  HAVE  ANY  PORTIONS  HAT  A^E  CONSTANT  TH  FT  A 
C  L^'NES. 

SMRF(  THETA  )  =  •  .  21  2 2 E7  07E  ^7377  TE’’ T  c  u ?- 5  2360: 
l«''  FORMAT  ("I  ••) 

I'T''  FORMATC  •* ,  "KGr  I D"  ,5  X ,  1  3 , 2  (5  X ,  ^1 E  .  1  C )  ) 

I  PORMAT  ("fGRI0",5X,13,2(5X,P10.1.  )  ) 

G'^VEN  THE  EOtJATICN  FOR  THE  POLE  SHOE  SURFACE,  CALCULATE  THE 
COOROINATES  FQP  THF  POINTS  A,  9,  C,  AVO  0. 

THE  CONSTRAINTS  AREI 

1)  THE  ARC  LENGTH  BETWEEN  POINTS  A  AND  9  ^UST  BE  60 
MECHANICAL  DEGREES. 

2T  THE  LINE  SEGMENTS  A3  AND  DC  MJST  PE  PARALLEL  TO  A 
LINE  THROUGH  THE  CENTER  OF  THE  GENERAruR  DIVIDING 
THE  POLE  IN  HALF, 

PI  =  3.  141Lg2ti>3F  8q7  93  23v3462b4  3  3S32H5D  0 

TAs7?;,'ior»pi/ieJ.^)D:' 

T0  =  1<=,  ’0  'PI/IEO.  TOO 
XY  =  OSORT(  j.EO.  *  3. 3bO:*3, 3000 
X7=0S0RT ( j.EOr  -  1, 3  20u^l.  320: ) 

XT  =  OSORT(C.*Di  ♦0.195  03*0.1950(5) 

XCC=yT+X7+XT 

YCC=X-^-XT 

RC  =  DSORT( (XCC'  XCC) ♦( YCC^YCC) ) 

R3  =  RC 

PA=SURF(TA> ♦R3 
RO=SURF(TO) ♦RP 

XA3=RA*0SIN(3i.a0  T*PI/18u.:0:) 

T3=('‘5.CDw»FI/13.i.CDv)+DASIN  (XA3/R3) 

'XC0=FD*0SIN  (3  -  .OO'j^P  1/180.00  i) 

TC=(A5.l,DC»  PI/13  r  .COu)-OASIN(XCl/RC) 

WRITE  (6,1 Vlt) 

PUNCH  THE  RGRIO  CARDS. 

RGRIO  points  are  TO  PE  EQUALLY  SPACED  ON  GIVEN  SEGMENTS, 

REMOVE  the  C  in  column  1  PRIOR  10  ALi.  HRITEIA,  ...  STATEMENTS 
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c  !«■  PUMCHEO  C«RO  OUTPUT  IS  OESIREO, 

Z  PUNCH  CftF'DS  FOR  SEGMENT  K»R  -  OMLT  THE  ?IGHTHAMO  ENT  POINT 

NG=LPG 

DO  1:  3=l,NPOMi 

AMG=  (TC/*  ’  (I-D'dej,  0  /PI) 

NG=NG+1 

WRIT'^  (m,1  1'.  )  )  G,RB,  ftNG 
0  WRITE  C^jl  u2(  )  T.G,R<3,ANG 

ir  CONTINUE 

3  OHMCH  '?GRTn  CAf^CS  FOR  SEGMENT  ^p  -  ON.Y  THE  RIGHTHANO  END  POINT 


N0H1=N0-1 

DO  2'.  I  =  1,NPM1 

ANG=  (TC+(  (  (TB-TC)  /NPMl)  ♦  (I-l)  ))  M  I3.j.  uDl'  /  PJ) 

NG=NG  *1 

WRITE  (3,i: Ir)  NG,P9, ANG 
:  WRITC- (7,1.  2' )  MG,RP,ANG 

?'  CONTINUE 

C  PUNCH  RGPin  CARDS  FOR  SEGMENT  RPA  -  BOTH  END  POINTS 

N0AN1=NPA-1 

00  3(  1  =  1, NPA 

AMG=(TB+(  ((  (PI/2.  '100) -T8)/NPAM1)*  (I-l )  > )  »  (  18  C.!  Or/PI) 

NG=NP-^1 

WRITE  (5,1 jIC)  NG,R9,ANG 
0  WRITE  (7, 1^,2' )  NG,RB,ANG 

71  CONTINUE 

Z  PUNCH  RGRIO  CARDS  FOR  SEGMENT  •'^3  -  ONLY  THE  RIGHTHAND  END  POINT 
C 

NG=LNG 

NPRH1=NF5-1 

00  1=1,NPBM1 

ANG= (TD/NFBMl) ' ( I-l) ♦ (18 0. (0 . /°I) 

NG=NG+1 

WRTTE(S,l(,ir)  NG,RO,ANG 
C  WRITE  (7,li 2^ )  NG,RO,ANG 

t”  CONTINUE 

0  PMNCH  RGRIO  cards  for  segment  K-  -  BOTH  END  POINTS 

NPM1=NF-1 

00  5.'  I  =  1,NFM1 

ANG=  (TO-K  ( (TA-TO)  /NPMl)  ♦  (I-l) ) )  •d  3  J.  ro:  /  PI) 

RAD=SUPF(ANG)*EB 
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f)  f)  f)  f*  Of''*’  f)  t>  O  O  f  1  t'  O  O  O  » 5  OC'it'i 


NG  =  NG  «•! 

WRITPCG,!  JIC)  ^'G,RAf),ANG 
WRITr(7,lCE:)  NG,RAD,ANG 
GONTINUr 


o’'VCH  CAh.rS  FOR  SEGMtNiT  \  -  ^DTH  E^'D  POINTS 


Ne’AMl  =NFf  -1 

on  6’  :=1,NFA 

(T  A+{  (  (  (PI/2.  .0.])-T  A)/M''AM1)  >■  ( I  -  1 )  )  )  M  18  .  T'’  /PT) 

NG=»4G  ♦I 

WRIT-  1  NG,PA,ANG 

WRITE  (■',  lv.2,  )  K'G,  RA,  AN'G 
<=  ■:  CONTI  MUE 


P"NCH  =;r,r.  to  CAFCS  FOR  SEGMENT  ><3 1  -  3DTH  END  POINTS 


NG1M1=NG1 -1 

DO  7.  I=1,NG1 

ANG=  (  {  (Pl/2.  .  0  )  /  NGIMD*  (I-l )  )•  (1  5’  .  .  D'  /  ®I ) 

ND  =  NG  +•! 

WRTTf  ('>,lcir  )  NG,RG1,ANG 
WRITE  (7,1.2"')  NG,RG1,ANG 
7Z  CONTINUE 

ALL  0«^  THE  RGRIO  CVROS  APE  NOW  PUNCHED 
PROCEED  TO  PUNCH  THE  TftIA  CARDS 

R'-^OV'"  rut  0  IN  COLUMN  1  PRIOR  TO  ALL  WRITE(T,  ...  STATEMENTS 
IF  PUNCH'^D  OUTPUT  IS  DESIRED. 

PUNCH  TRIA  CARDS  FCR  REGION  T1  (DETWEEN  THE  DUADRATURF  AVIS 
AND  THE  RIGHT  SIDE  OF  THE  POLE  SHOE 

KEY=1 
NT=LNT 
KT=1 
<D  =  1 

DNFPP=pL0AT (NFP) /FLOAT(NPB) 

P"  0'<TD=FLCAT  (KT) /FLOAT  (KP) 

IF  (nKTB.  GT  .PNFpq)  GO  TD  9t' 

IF  (kt.eo.nfb)  go  TO  iro 
CALL  TRIAl  (RT,K3,LNG,LPG,Nr jKEY) 

GO  TO  PC 

"P  CALL  TRIA2  (KT  , K B , LNG ,L P6,NT , KEY ) 

GO  TO  S'* 

PUNCH  TRIA  CAROS  FOR  REGION  T2.  (THE  POLE  SHOE) 
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Cl  CJ  O 


1  ' 

LNG=L  NG+NFF-1 
L0G=LPG+NP9-1 
<T  =  1 
<?  =  ! 

n  JFP=PLC4T(rF) /FLOAT  (NP) 
l:  '  KTl^PLOAT  (KT) /c^LOAT  (KD) 

C^x'TP.GT.ONrs)  GO  TO  12 

C^T.EO.LF)  go  to  13' 

TILL  TPIfl  (<T,KF,lNG,LPG,KT, KEY) 

GO  TO  il‘ 

1?'  call  TP1A2  (<T,<e,LNG,LPG,wT,Kt:Y) 

GO  ^0  11- 

oiJNCH  TRIA  CAFOS  "EGION  T3  (OETMEEN  THE  OUADFATJPt  AXIS 

AMO  THE  left  SIPL  CF  TH£  POLE  SHOE) 


17"  i<EY=t 

LNG=LNG+NF-1 
L3G=LPG  +  ‘  P-1 
Kr  =  i 
<”1  =  1 

ONFPAsFLOAT (NFA) /FLOAT(NPA) 

I'-’  0<T3=FLOfiT(t<T)/FLCAT(k'P) 

IP  (DKTE.GT  .0^'FPA)  GO  TO  13J 

IF  (KT.EO.NFC)  GO  TO  If*, 

CALL  TPltl  (KT,K9,LNG,LPG,Nr,Kr:Y) 

GO  TO  lAL 

IT  CALL  TRIA2  (KTjKOjLNGjLPGjHTjKEY) 

GO  TO  IT 

P'INCH  T^IA  Cards  for  region  T4,  (lower  half  of  the  CIR  GAO) 

IF*'  KEY=1 

LNG=LNG*NFA 
LoG=LNG-NFA-NF-NF5+2 
KT  =  1 
<0  =  1 

DNGFM2  =  FL0AT  O.Gl)  /FL  OAT  (  NFA+MF4-Hr  R-2) 

177  d<T3=FL0AT (KT) /FLOAT (K9) 

IF  (OKTP,  GT.nNCF^2)  GO  TO  13.' 

IF  (KT.EO.NGl)  GO  TO  193 

CALL  TRIAl  (KT,<0,LNG,LPG,NT,  KEY) 

GO  TO  in 

133  CALL  T0IA2  ( KT , < 9 ,LN G ,L PG ,NT , <EY) 

GO  TO  171 

PUNCH  TRIA  CARDS  FOR  REGION  T5,  (UPPER  HALF  OF  THE  AIR  GAP) 
IPn  KEY=1 
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r 


LN';=LMS 

L°G=LPG*(NFA+NF+NFB-2) 

KT  =  l 
<^=1 

(N£)  /FLOAT  (NGl) 

?  ■  n'<TT  =  '^LCST  (XT) /float  (<R) 

I*"  (ni<TB,GT.DK‘'Gl)  GOTO  21  ■ 

I-  {KT.EO.^.’S)  GO  TO  22: 

CALu  T^IAl  ,<'’,LN6,LPG,MT,  <cT> 
GO  TO  2vr 

?' •  CALL  •^=?IA2  (KT  ,><?  ,LNG,LPG,^'^  jKTY) 
GO  TO  2i  L 
2?^  STOP 
ENO 
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J 

Z  SUBt?OUTINE  PUNCHES  TRIA  CARJS  ED?  TrPE  1  EIMTE  ELEMENTS. 

C 

C  S“-  TH'STSl  FO'  AN  EXPLANATION  OF  £l«  \/ARIARLES. 

<./ 

C  -■‘^OVT  the  C  IV  CCLUMN  1  PrlC-J  TJ  Al.  WRiTrK?,  .  .  .E  T  AT  EH  "NTS 
C  !'■  O!JNCH'^0  OUTPUT  IS  OESlRfO. 

•J 

Z  Tver  ^  PTNITr  ELEI-^NTS  AkE  TRIANGLES  WITH  OOWNWAPP  POINTING 
C  fiP-XES. 

S'IPROUTIKC  TRIAI  (  KT  jKPfLNGjLPGjNr  ,  key  > 

NPrK" +LNG 
N1  =  N  ?  +  l 
NS=<S  +LPG 
KT=KT ♦I 
NTrsl”  +1 

IF  (KEv.EO.l)  GO  TO  2'' 

W^IT-  (F:,l  .  1'  )  NT,N1,N2,N3 
:  WRIT-  (7,1  ,11)  VT,N1,N2,N3 

I’l*  FORiiTC*  "j**?"  1A",-+(EX,  13)  ) 

l''ll  PORMAT  (••TKIA",L  (:=X,I3)) 

RETUP  N 

J 

2'  WRITE  (f-,,1',2:  )  NT,N1,N2,N3 

C  WRITE  (7,i:21)  NT,t.l,N2,N? 

1^2)  FORMATC*  ••,"TFIA*‘,k(&X,I3),5X,"AIR") 
i'’21  FOR.MA‘^("TKIA”,t(rx,l3)  ,5X,"AIP") 

PET'JPN 

ENT 
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T'JIS  SU^^OUTINE  PUNCHES  TPIA  CftR^S  FC?  TlfPE  ?  FJNIT-  ELEMENTS 


C  E'C  THESISl  FO=  AN  EXFlANATION  OF  A^L  ^A^IA^.ES. 

0  ^  ‘'T/-  T>^!r  ^  TK'  CCLUP.N  1  PRIOri  Tn  A;..  . . .  ^  t  t ;  th  r^JT5 

:  ^'.inCHFO  OUTPUT  IS  DESIRED. 


0  -Y'j-  FTNI'^E  cLcN^NTS  ARE  TRIANGLES  WITH  L'PWA.-C  Pf^’M^ING 
C  Ac-x;.<;, 

S'l'^ROUTINt  TR1A2  (KTjKBjLNGjLPGjNT,  KEY) 

Nl=KT ♦LNG 
N2=<p+LPG 
N'^=S2  +  1 
<q  =  <R  4-1 
NT=NT  +1 

ir  (KEY.f^.l)  GO  TO  2:? 

WRITE(S,1.H)  MT,Nl,N2,N3 
0  WRITE  ('’,l..ll)  NT,N1,(m2,N3 

IM'  FORMATC* 

imi  FORMAT<”TRiA",A  (5X,I3)) 
return 

J 

2''  WRIT-' (G,l -2t  )  NT,N1,N2,N3 
:  WRirE(7,l. 21)  NT,N1,N2,N3 

I’O''  FOPIA'^C*  ••,"THA",R(yX,I3),5X,"AIF”) 

1''21  FORMAT  ("TPl  A",l-  (f^X,I3)  ,5X,"AIR") 

RET'JPN 

ENO 
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Cecil  Stucrko  was  born  2h  September  1949  in  Sedalia, 
Missouri.  He  grew  up  in  Kansas  City,  Kansas  where  he 
graduated  from  Wynadotte  High  School  in  1967.  He  attended 
Kansas  City,  Kansas  C.’ommunity  Junior  College  tiie  following 
year.  He  tlien  attended  Kansas  University  from  1968  to  1971. 
Short  of  graduation,  lie  enlisted  in  the  United  States  Air 
force.  After  serving  one  year  at  Richards  Ccrhaur  Air  force 
base,  Grandview,  Missouri,  lie  returned  to  Kansas  University 
under  the  Airman  liducation  and  Commissioning  Program.  He 
received  Bachelor  of  Science  Uegrees  in  lilcctrical  lingincering 
and  lingincering  Pliysics  from  Kansas  University  in  1978.  He 
attended  Officer's  Training  School  and  received  an  Air  I'orcc 
Commission  on  24  March  197.S.  He  was  tlien  assigned  to  the 
Air  force  Ac ro - Projui  1  s  ion  Labortatory  at  Wright -Patterson 
Air  force  Base  until  entering  the  School  of  lingincering.  Air 
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hcL-ausc  the  rotor  MMl-  lias  only  the  velocity  of  the  rotor. 

The  MMlh^  j-ji-c.^cnt  at  all  tlio  .stator  slots  of  an  arma.ture  coil 
group  are  reflected  into  a  single  slot  in  an  attempt  to  control 
the  air  gap  flu.\  density  harmonics  he  selectively  positioning 
increments  of  rotor  reluctance.  This  method  fails  because  a  gen¬ 
eral  solution  for  tiie  reflected  flu.x  densities  could  not  he 
o  h  t  a  i  n  e  d  . 

Sj^ecifying  the  air  gap  reluctance  incrementally  a.s  tlic  rotor 
moves  under  successive  stator  slots  is  also  tried  as  a  means  of 
harmonic  control.  Tliis  method  proves  unworkable  due  to  negative 
.\IMF  reciui  rement  s  and  conflicting  geometry  specifications. 

Finite  FMcmcnt  analysis  is  introduced  as  a  computer  aid  to 
the  design  of  a  magnetic  structure. 

Recommendations  for  further  studv  are  included. 
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